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ADAPTIVE  BEAMFORMING  ALGORITHMS  FOR 
HIGH  RESOLUTION  MICROWAVE  IMAGING 

Bernard  D.  Steinberg 
Bongsoon  Kang 


0 .  INTRODUCTION 

This  report  is  an  algorithmic  study  of  adaptive  beam¬ 
forming,  which  is  a  means  for  self-calibrating  a  radio  cam¬ 
era.  A  radio  camera  is  a  bistatic  microwave  radar  with  a  huge 
antenna  array,  the  array  being  sufficiently  large  to  enable 
the  system  to  achieve  the  fine  angular  or  cross-range  resolu¬ 
tion  that  hitherto  had  been  available  only  with  optical  or 
near-optical  instruments. 

The  bulk  of  the  work  performed  to  date  falls  into  two 
categories.  In  the  first  body  of  work  the  self-calibrating 
portion  of  the  phased  array  system  is  driven  by  the  reradia¬ 
tion  from  a  prominent  reflector  having  large  radar  cross  sec¬ 
tion  and  small  physical  size.  The  conditions  under  which  such 
a  phase-synchronizing  source  operates  satisfactorily  are 
given  in  [8],  [4].  Algorithms  of  the  second  class  are 
designed  to  operate  upon  the  backscatter  from  statistically 
homogeneous  targets  and  terrain  [18],  [17].  Examples  of  such 
backscattering  sources  are  sea  clutter,  crowns  of  trees  in  a 
forest  and  large  expanses  of  sand. 
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Because  of  the  lessened  dependence  upon  specificity  of 
source  distribution  the  algorithms  in  the  second  body  of  work 
are  generally  referred  to  as  scene-free  algorithms.  The  RADC 
support  for  the  reported  work  focused  upon  the  scene-free  ABF 
procedures . 

For  completeness,  however,  both  classes  of  algorithms 
are  described.  Section  1  is  introductory,  describing  uses  of 
the  radio  camera . Section  2  gives  a  brief,  partial  history  of 
self-calibration,  discusses  certain  interesting  properties  of 
ABF  algorithms  and  indicates  desirable  properties  of  such 
algorithms.  The  two  classes  of  algorithms  are  described  in 
Sections  3  and  41.  These  sections  also  show  experimental 
results  of  the  application  of  these  algorithms  to  high  reso¬ 
lution  X-band  (3  cm  wavelength)  radar  data  obtained  from  VFRC 
radio  cameras.  Section  5  compares  the  relative  computation 
times  required  by  algorithms  from  both  classes.  Section  6  is 
the  Summary.  The  Appendix,  authored  by  Ph.D.  candidate 
Bongsoon  Kang,  presents  additional  experiments,  analytical 
details  and  further  comparisons  between  algorithms. 

It  is  found  experimentally  that  highly  distorted  phased 
arrays  can  be  self-calibrated  under  a  very  wide  range  of  con¬ 
ditions.  The  dominant  scatterer  algorithm  (DSA) ,  which  is 
the  basic  procedure  in  the  first  class,  almost  always  can  be 
made  to  self-calibrate  on  man-made  structures.  A  modifica- 

^hese  sections  form  the  bases  for  two  chapters  in  a  forthcoming  book  to 
be  published  by  John  Wiley:  B.  D.  Steinberg  and  H.  Subbaram, 

Microwave  Imaging  Techniques. 
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tion,  the  multiple  scatterer  algorithm  (MSA) ,  extends  its 
usefulness  under  conditions  of  poor  target  to  clutter  ratio. 
The  scene-free  ABF  algorithms  are  superior  when  the  source 
distribution  is  statistically  homogeneous,  e.g.,  sea  clutter. 
Although  the  scene-free  techniques  were  designed  primarily 
for  statistically  homogeneous  backscattering,  the  experimen¬ 
tal  evidence  indicates  that  they  perform  exceedingly  well 
under  those  conditions  in  which  the  dominant  scatterer  algo¬ 
rithms  also  work.  Thus  the  scene-free  algorithms  have  been 
found  to  be  more  generally  applicable  than  the  dominant  scat¬ 
terer  algorithms,  they  do,  however,  require  greater  computa¬ 
tional  strength.  However,  because  of  the  continuing  improve¬ 
ment  over  the  years  in  speed,  capacity  and  cost  of  realtime 
signal  processing  devices,  it  is  concluded  that  the  design  of 
future  radar  systems  requiring  very  large  antennas  should  be 
based  upon  the  scene-free  algorithms. 

The  authors  are  grateful  to  Dr.  Robert  Shore  of  RADC, 
Hanscom  AFB,  Bedford  MA  for  his  careful  editing  of  the  draft 
of  this  report. 
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1.  THE  RADIO  CAMERA 

The  diffraction-limited  resolution  of  an  aperture  (lens,  te¬ 
lescope  mirror,  microwave  dish,  phased  array,  ultrasonic  trans¬ 
ducer)  is,  in  radians,  approximately  the  reciprocal  of  the  size  of 
the  aperture  in  units  of  wavelength.  When  the  aperture  is  too 
large  to  maintain  the  rigidity  necessary  for  diffraction-limited 
performance  (approximately  one-tenth  of  a  wavelength)  it  is  neces¬ 
sary  to  introduce  a  low-order  intelligence  into  the  antenna  system 
so  as  to  dynamically  self-calibrate  it  as  it  undergoes  distortion. 

Such  an  aperture,  i.e.,  one  too  large  to  maintain  its  shape, 
is  called  a  large  aperture.  Large  apertures  include  distributed 
antennas,  that  is,  antennas  the  parts  of  which  are  located  on  sep¬ 
arate  and  independent  platforms.  Examples  of  such  apertures  and 
their  potential  applications  are  shown  in  Figures  1.1  to  1.4. 

Figure  1.1  depicts  either  an  air  traffic  control  (ATC)  envi¬ 
ronment  or  an  aircraft  landing  problem.  An  airport  is  shown  with 
an  ATC  radar  illuminating  a  target  with  a  microwave  pulse  and 
receiving  an  echo.  The  echo  is  displayed  as  a  blip  on  the  monitor 
scope  in  range-angle  coordinates.  The  presence  of  the  blip  consti¬ 
tutes  detection  and  the  coordinates  indicate  the  distance  to  the 
target  and  its  bearing.  No  shape  information  is  obtained  because 
the  width  of  the  beam  is  very  large  compared  to  the  size  of  the 
airplane,  as  is  easy  to  see  from  the  following  calculation. 
Diffraction  theory  teaches  that  the  angular  width  of  the  beam,  ir. 
radians,  is  approximately  the  ratio  of  the  wavelength  of  the  radi- 
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FIGURE  1.1.  AIRPORT  APPLICATION  OF  THE  RADIO  CAMERA.  A  CONVEN¬ 
TIONAL  RADAR  TRANSMITTER  ILLUMINATES  THE  AIR¬ 
PLANE  WITH  A  PULSE.  ECHOES  ARE  RECEIVED  BY  THE  RADIO 
CAMERA  MODULES.  EACH  OF  WHICH  CONSISTS  OF  AN  AN¬ 
TENNA  AND  A  COHERENT  MICROWAVE  RECEIVER.  THE  MOD¬ 
ULES  ARE  DISTRIBUTED  ACROSS  THE  AIRPORT  WHEREVER 
THERE  IS  EMPTY  SPACE  THE  ECHOES  ARE  TRANSMITTED  TO 
THE  RADIO  CAMERA  SIGNAL  PROCESSOR  WHICH  CREATES  A 
HIGH  RESOLUTION  IMAGE  OF  THE  TARGET.  THE  BEAM- 
WIDTH  OF  THE  RADIO  CAMERA  IS  lO^irad  WHEN  THE  WAVE¬ 
LENGTH  IS  5cm  AND  THE  AIRPORT  IS  5km.  (FROM  [4].) 
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ation  to  the  size  of  the  antenna.  A  typical  value  of  this  ratio  is 
1/30.  The  beam  cross-section  at  a  target  distance  of,  say,  30  km 
is  one  km  whereas  the  size  of  a  large  airplane  is  only  40  to  50  m. 
Because  the  target  is  so  much  smaller  than  the  beam  cross-section 
it  reflects  as  a  single  point  and  appears  as  such  in  the  radar 
display . 

The  airport,  on  the  other  hand,  is  several  km  in  size. 
Imagine  that  an  imaging  system  consisting  of  hundreds  of  small 
antennas  and  receivers,  plus  the  proper  central  signal  processor, 
is  spread  out  over  the  airport  in  the  manner  shown  in  Figure  1.1. 
Assuming  that  the  signal  processor  could  appropriately  process  the 
echoes,  the  width  of  the  beam  that  would  be  formed  would  be  deter¬ 
mined  by  the  size  of  the  airport  rather  than  the  size  of  each 
individual  antenna.  For  example,  let  the  airport  be  5  km  and  the 
wavelength  5  cm.  The  beamwidth  of  the  system  would  be  1/100,000 
rad  and  the  beam  cross-section  at  the  target  would  be  30  cm.  This 
30  cm-cross-section  beam  could  then  be  scanned  across  th.  _arget 
to  provide  an  image  with  more  than  100  resolution  cells  in  the 
cross-range  or  azimuth  direction.  Experience  with  very  high  reso¬ 
lution  imaging  radar  at  the  Valley  Forge  Research  Center  (VFRC)  of 
the  University  of  Pennsylvania  shows  that  at  most  30  to  50  resolu¬ 
tion  cells  on  a  target  are  adequate  to  define  the  target.  Hence  in 
this  ATC  example  the  radar  controller  would  be  provided  with  ade¬ 
quate  resolution  in  the  cross-range  direction  to  differentiate 
between  types  of  aircraft. 

Of  course  high  resolution  in  the  range  direction  is  also  nec¬ 
essary  for  it  provides  the  second  dimension  needed  for  two-dimen- 
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sional  imaging.  High  resolution  in  range  is  obtained  through  the 
use  of  a  wide  bandwidth  transmitted  signal.  In  the  VFRC  experi¬ 
ments  shown  later  in  this  report  a  7  ns  pulse  having  150  MHz  band¬ 
width  is  employed.  With  such  a  waveform  a  range  resolution  of  1  m 
is  achieved. 

An  imaging  instrument  of  this  kind  is  called  a  radio  camera 
[1],  It  is  basically  a  bistatic  imaging  device  operating  at 
microwave  wavelengths.  Like  a  conventional  optical  camera  it  oper¬ 
ates  upon  the  reflections  from  an  illuminator  that  does  not  have 
to  be  colocated  with  the  receiving  aperture.  There  the  similarity 
stops,  however.  Not  only  must  the  receiving  aperture  be  large 
enough  to  provide  the  desired  resolution,  but  because  the  wave¬ 
length  of  the  radiation  is  three  to  six  orders  of  magnitude  larger 
than  optical  radiation,  there  are  several  fundamental,  inherent 
design  differences. 

These  differences  fall  into  three  categories.  The  first  is 
aperture  design.  The  second  category  is  the  scattering  properties 
of  targets  vs.  wavelength.  The  third  is  propagation  and  weather 
effects . 

The  first  category  is  the  most  critical.  Because  of  the  huge 
aperture  size  required  at  microwaves  to  achieve  the  modest  angular 
resolution  of  common  optical  instruments  the  radio  camera  aperture 
must  be  assumed  by  the  designer  to  be  intrinsically  distorted. 
That  is,  the  designer  is  forced  to  assume  that,  no  matter  how 
carefully  s/he  specified  the  antenna  element  locations  and  their 
complex  weights,  the  actual  antenna  will  differ  sufficiently  so  as 
to  destroy  the  desired  radiation  pattern. 
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Consider  the  illustration  above.  The  5  km  radio  camera  oper¬ 
ating  at  5  cm  wavelength  is  as  large  in  units  of  wavelength  as  a 
50  mm  optical  telescope.  While  the  latter  can  be  constructed  from 
a  small  lens  purchased  from  a  scientific  supply  shop,  the  former 
cannot  be  designed  from  standard  antenna  design  practice  alone. 
First,  it  cannot  be  assumed  that  the  required  positional  tolerance 
of  less  that  1/10  wavelength  can  be  achieved  when  the  antenna  ele¬ 
ments  are  deployed  over  100,000  wavelengths.  Second,  when  the 
receivers  are  spaced  so  far  apart  their  electromagnetic  couplings 
to  their  local  environments  will  differ.  The  consequence  is  a  ran¬ 
dom  component  to  the  driving  point  impedance  of  each  antenna  ele¬ 
ment  in  the  giant  array,  which  appears  as  a  random  complex  gain  or 
weight  in  each  signal  channel  of  the  imaging  system.  Third,  the 
system  will  be  paced  by  a  local  oscillator  distribution  subsystem 
that  must  remain  frequency-  and  phase-locked.  With  a  system  of 
radio  camera  size  it  must  be  assumed  that  oscillators  will  drift 
and,  similarly,  that  electronic  circuits  in  the  receivers  will 
change  their  complex  gains  with  age.  Fourth,  impedance  mismatches 
throughout  the  system  will  introduce  further  random  amplitude  and 
phase  errors  in  the  signal  channels.  Fifth,  these  random  random 
errors  can,  and  in  all  likelihood  will,  vary  with  time. 

In  principle  each  of  these  distorting  effects  can  be  tuned 
out.  In  practise  dynamic  self-calibration  is  required.  Such  tech¬ 
niques  have  been  developed  at  VFRC  during  the  last  several  years 
and  the  comparison  between  the  performances  of  several  of  them  is 
the  subject  of  this  report. 
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The  second  category  of  difference  between  radio  camera  design 
and  the  design  of  optical  instruments  is  due  to  the  vast  differ¬ 
ences  in  the  scattering  properties  in  the  two  wave  bands.  Most 
manmade  surfaces  are  rough  to  optical  radiation  and  therefore 
broadly  scatter  optical  energy  in  all  those  directions  not  prohib¬ 
ited  by  geometric  shadowing.  By  and  large  those  same  surfaces  are 
smooth  to  microwave  radiation  and  reflect  more  or  less  specularly. 
Thus  most  surfaces  appear  black  in  a  microwave  image,  and  when 
they  do  show  up  they  do  so  as  highlights  or  specularit ies .  As  a 
consequence  the  dynamic  range  is  50  to  90  dB  in  a  microwave  image, 
whereas  it  is  only  20  to  30  dB  in  an  optical  scene  even  in  bright 
sunlight  on  a  beach  or  in  snow.  Weak  scattering  from  edges  and 
corners  of  targets,  which  tend  to  define  their  shapes,  become  ob¬ 
scured  by  thCT  strong  specular  returns. 

Another  problem  resulting  from  the  prevalence  of  specular 
reflections  is  that  their  sidelobe  responses  coherently  interfere 
to  produce  a  random  speckle  pattern  in  the  image.  The  effect  is  an 
unrealistic  "texture"  to  the  target  image.  Still  another  problem 
resulting  from  destructive  interference  is  a  random  cancellation 
of  image  signals  at  the  edges  or  peripheries  of  targets.  The 
effect  is  to  make  large  targets  appear  as  groups  of  smaller  ones. 

Several  techniques  are  applicable  here.  One  is  diversity  com¬ 
bining  of  images  obtained  from  radiation  patterns  with  statisti¬ 
cally  independent  side  radiation  patterns.  A  second  is  image 
deconvolution  to  remove  the  specular  returns.  A  third  is  post- 
image-formation  digital  pixel  processing.  This  topic  is  not  cen¬ 
tral  to  this  report,  however,  and  will  not  be  discussed  further. 
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The  third  category  of  difference  in  design  between  microwave 
and  optical  imaging  systems  is  in  propagation  and  weather  effects. 
The  atmosphere  is  not  wholly  transparent  in  either  waveband. 
Internal  waves  in  the  atmosphere  induce  local  variations  in  the 
refractive  index,  which  distort  wavefronts.  For  an  aperture  to 
perform  at  its  diffraction  limit  its  size  must  be  sufficiently 
small  for  a  wavefront  arriving  at  it  to  appear  planar.  That  is, 
its  size  must  be  less  than  the  correlation  distance  of  the  refrac¬ 
tive  index  variations  of  the  propagation  medium  at  the  wavelength 
at  which  it  operates.  The  lens  size  in  comm:  optical  instruments 
such  as  the  eye,  the  camera  and  the  small  telescope  meet  this  con¬ 
dition,  although  larger  telescopes  such  as  the  200  in  Hale  tele¬ 
scope  on  Mount  Palomar  do  not.  It  is  questionable  whether  a  5  km 
microwave  antenna  would  satisfy  this  requirement.  Failure  to  meet 
this  condition  means  that  the  antenna  array  is  confronted  with  a 
phase-perturbed  wavefront  and  the  measured  signals  in  each 
receiver  channel  has  some  unknown  phase  error  caused  by  it.  A 
dynamic  self-calibrator  of  the  array  would  also  solve  this  prob¬ 
lem  . 

In  addition  to  propagation-induced  phase  errors  an  array  as 
large  as  the  one  under  discussion  is  likely  to  experience  differ¬ 
ent  weather  conditions  across  it.  The  nominal  size  of  a  rain  cell 
is  about  5  km.  Consequently,  one  portion  of  the  array  could  be  in 
sunlight  while  it  could  be  raining  in  another  part  of  the  array. 
Water  pouring  over  the  radomes  covering  some  of  the  receivers 
could  induce  randorti  phase  shifts  to  the  waves  passing  through. 
Again,  dynamic  self-calibration  solves  this  problem. 


Thus  the  radio  camera  is  a  bistatic,  high  resolution  imaging 
system  with  a  large,  intrinsically  distorted  antenna  and  electri¬ 
cal  system,  both  of  which  are  dynamically  self-calibrated .  The 
generic  term  for  the  various  self-calibration  algorithms  is  adap¬ 
tive  beamforming  (ABF) . 

Figures  1.2  to  1.4  illustrate  other  uses.  Figure  1.2  shows 
an  airborne  radio  camera.  The  antenna  elements  are  flush-mounted 
in  the  skin.  Notwithstanding  the  vibration  of  the  skin  and  the 
flexure  of  the  airframe  the  ABF  subsystem  in  the  image  processor 
compensates  for  the  distortion  and  the  output  image  is  from  a 
diffraction-limited  system.  The  interesting  feature  of  the  air¬ 
borne  radio  camera  is  that  large  jet  aircraft  are  as  large  as  the 
human  eye  when  measured  in  units  of  wavelength.  Consequently  the 
image  displayed  to  the  navigator  would  have  the  same  angular  reso¬ 
lution  as  the  human  optical  system.  Such  a  passing-scene  display 
with  human  optical  resolving  power  would  be  available  to  the  air¬ 
man  at  night  as  well  as  during  the  day  and  under  all  weather  con¬ 
ditions  . 

Figure  1.3  depicts  a  battlefield  surveillance  problem.  A  dis¬ 
tant  transmitter  on  a  hill  is  shown  illuminating  a  target  region 
of  interest.  A  large  receiving  array  is  hurriedly  deployed,  per¬ 
haps  from  the  tailgate  of  a  truck.  Phase  and  frequency  synchro¬ 
nization  information  is  radioed  from  the  transmitter  to  each  of 
the  receiver  modules.  The  locations  of  the  free-standing  microwave 
receivers  are  known  only  crudely.  Yet  ABF  forms  and  focuses  a  beam 
in  the  target  region  and  scans  the  beam  to  obtain  a  high  resolu¬ 
tion  image.  If  the  target  region  is  10  km  from  the  receiving  array 
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FIGURE  1.2.  AIRBORNE  APPLICATION  OF  THE  RADIO  CAMERA.  ANTENNA  ELEMENTS 
ARE  DISTRIBUTED  THROUGHOUT  THE  AIRFRAME  TO  FORM  A  NONRIGID 
PHASED  ARRAY.  THE  ELEMENTS  RECEIVE  GROUND  ECHOES  FROM  THE 
TRANSMITTER  IN  THE  NOSE  RADAR.  THE  SIGNAL  PROCESSOR  COHERENTLY 
COMBINES  THESE  ECHOES  TO  FORM  A  HIGH  RESOLUTION  IMAGE  OF  THE 
PASSING  SCENE  BELOW  THE  AIRPLANE.  BECAUSE  THE  AIRPLANE  IS 
APPROXIMATELY  AS  LARGE  AS  THE  HUMAN  EYE  (AS  MEASURED  IN 
WAVELENGTHS)  THE  RESOLUTION  IN  THE  PASSING  SCENE  IS  COMPARABLE 
TO  HUMAN  VISION  AND  IS  AVAILABLE  NIGHT  AND  DAY  AND  IN  ALL 
WEATHER.  (FROM  [4].) 
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FIGURE  1.3.  A  BATTLEFIELD  SURVEILLANCE  APPLICATION  OF  THE  RADIO 
CAMERA  ILLUSTRATING  ITS  BISTATIC  NATURE.  THE  TRANS¬ 
MITTER  IS  SEPARATED  FROM  THE  RECEIVING  ARRAY  AND 
ILLUMINATES  THE  TARGET  AREA  TO  BE  IMAGED.  THE  ARRAY 
IS  HASTILY  DEPLOYED  AND  THEREFORE  IS  ILL-SURVEYED 
AND  BADLY  DISTORTED.  THE  RADIO  CAMERA  FORMS  A  DIF- 
FRACTION-LIMITED.  HIGH  RESOLUTION  IMAGE  OF  THE 
SCENE  (FROM  [4].) 
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and  it  is  desired  to  obtain  aim  cross-range  resolution,  the 
array  must  be  10,000  times  the  wavelength.  Assuming  X-band  (3  cm) 
transmission  the  modules  need  be  deployed  over  300  m. 

There  is  a  highly  favorable  theoretical  relationship  between 
the  size  of  the  angular  sector  (called  the  field  of  view  (FOV)  ) 
that  may  be  scanned  and  the  initial  error  or  uncertainty  in  the 
locations  of  the  antenna  elements.  Let  Gbe  the  standard  deviation 
of  element  position  error  and  let  A.  be  the  wavelength.  Then  the 

allowed  position-error  tolerance  is  [9] 
a  =  A./2tc(FOV) 

With  conventional  design  the  tolerance  is  approximately  A./4 n. 
Therefore  the  tolerance  in  element  placement  is  eased  or  loosened 
by  the  factor  2/FOV.  In  general,  the  FOV  need  be  no  larger  than 
the  target  being  imaged.  In  the  ATC  example  of  Figure  1.1  the  min¬ 
imum  FOV  =  30  m/3  km  =1/1000.  This  means  that  the  tolerances  (in 
manufacturing,  installation,  operation  and  maintainence)  are  eased 
by  nearly  three  orders  of  magnitude  in  that  example . 

It  is  evident  from  the  example  shown  in  Figure  1.3  that  the 
individual  receivers  of  the  radio  camera  could  be  on  separate 
platforms.  The  platforms  could  be  ground  based,  airborne  or  orbit¬ 
ing  stations  in  space.  Figure  1.4  illustrates  a  space-borne  radio 
camera  with  loosely  connected  or  tethered  receivers.  The  central 
vehicle  carries  the  transmitter  in  a  low  orbit  which,  for  the  sake 
of  illustration,  is  taken  to  be  at  a  height  of  200  km.  After 
unfurling  from  the  station  the  receiver-carrying  arms  form  a 
highly  nonrigid  10  km  array.  If  the  wavelength  is  10  cm  the  array 
is  100,000  wavelengths  and  the  beam  cross-section  on  earth  is  2  m. 
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FIGURE  1.4.  SPACE-BASED  RADAR  SURVEILLANCE  APPLICATION  OF  THE 
RADIO  CAMERA.  A  RADAR  TRANSMITTER  IN  THE  SATELLITE 
ILLUMINATES  THE  GROUND.  RECEIVERS  IN  ARMS  UNFURLED 
FROM  THE  SATELLITE  DETECT  ECHOES  FROM  EARTH.  NOT¬ 
WITHSTANDING  THE  DYNAMIC  SHAPE  VARIATION  OF  THE 
RECEIVING  ARRAY  THE  RADIO  CAMERA  SIGNAL  PROCESSOR 
FORMS  A  COHERENT,  DIFFRACTION-LIMITED  IMAGE  OF  THE 
AIR  TRAFFIC  BELOW.  GIVEN  THE  ORBIT  ALTITUDE  AS  200 
km,  THE  ARRAY  SIZE  AS  10  km  AND  A  10-cm  WAVELENGTH, 
THE  BEAMWIDTH  OF  THE  RECEIVING  ARRAY  IS  10  urad  AND 
THE  BEAM  CROSS-SECTION  ON  THE  EARTH  IS  2  m. 
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This  remarkable  resolving  power  is  potentially  available  from  a 
space-borne  radio  camera. 
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2.  ADAPTIVE  BEAMFORMING  (ABF) 

Of  the  many  problems  encountered  in  long  distance  microwave 
imaging,  one  problem,  that  of  aperture  distortion,  is  in  a  cate¬ 
gory  by  itself:  imaging  with  large  apertures  is  not  possible 
unless  the  problem  is  solved.  As  discussed  in  Section  1,  the  dis¬ 
tortion  is  both  geometric  and  electrical.  These  distortions  are 
intrinsic  to  all  systems  sufficiently  large  to  provide  the  resolu¬ 
tion  required  of  long  distance,  long-wavelength  imaging  systems. 

Two  sets  of  self-calibration  procedures  are  described  in  the 
following  two  sections.  Adaptive  beamforming  (ABF)  is  the  generic 
term.  Other  names  are  self-cohering,  self-calibrating,  self-focus¬ 
ing,  auto-focusing,  self-phasing  and  self-synchronizing. 

Adaptive  control  of  the  weight  vector  of  a  large  phased  array 
is  in  experimental  use  today  in  terrestrial  high  resolution 
microwave  imaging  as  well  as  in  radio  astronomy.  The  radio  camera 
was  first  described  in  a  concept  paper  in  1973  [1]  .  The  first 
radio  camera  experiment  demonstrating  ABF  was  performed  two  years 
later  [2]  .  The  algorithm  upon  which  most  subsequent  radio  camera 
experiments  has  been  based  was  published  in  1981  [3],  with  further 
details  in  1983  [4]  .  The  work  reported  in  these  and  other  refer¬ 
ences  demonstrate  that  diffraction-limited  imagery  can  be  obtained 
from  very  large  antenna  arrays  even  when  they  are  highly  dis¬ 
torted  . 

Since  the  early  part  of  this  decade  the  same  principle  has  been 
used  in  radio  astronomy  to  correct  for  uncertainties  in  the  loca- 
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tions  of  the  phase  centers  of  the  dishes  that  form  large  radio 
astronomy  arrays,  as  well  as  to  compensate  for  wavefront  distor¬ 
tion  induced  by  spatial  variations  in  the  refractive  index  of  the 
propagation  medium.  Review  articles  by  Cornwell  for  the  use  of 
these  procedures  in  radio  astronomy  are  found  in  [5]  and  [6] . 

An  interesting  property  of  the  radio  camera  is  that  it  is  a 
near-field  instrument.  In  actuality  its  means  of  operation  and  its 
performance  are  independent  of  whether  the  target  is  in  the  near 
field  or  the  far  field.  In  practice,  however,  high  resolution 
imaging  implies  near-field  operation.  The  transition  between  the 
near  field  and  the  far  field  occurs  at  that  distance  at  which  the 
wavefront  from  a  source,  when  viewed  at  the  antenna  array,  changes 
from  a  quadratically  curved  front  to  an  approximately  planar 
front.  This  distance  is  in  the  neighborhood  of  L2A,  where  L  is  the 

size  of  the  aperture.  At  this  distance  the  cross-range  dimension 
of  the  resolution  cell  equals  L.  At  greater  distances  the  cross¬ 
range  dimension  is  proportionally  larger  than  L  and  at  shorter 
distances  it  is  proportionally  smaller,  assuming  that  the  array  is 
properly  focused. 

An  easy  test  of  near-  vs.  far-field  operation  of  the  radio  cam¬ 
era  is  whether  or  not  the  cross-range  resolution  is  smaller  than 
or  exceeds  the  antenna  size.  It  is  easy  to  see  that  for  practical 
radar  usage  the  imaging  system  must  automatically  and  dynamically 
accommodate  to  near-field  operation.  The  minimum  available 
beamwidth  is  \/L  and  the  minimum  available  cross-range  resolution 

for  a  system  focused  to  range  R  is 
AS  -  A.R/L 
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A  typical  value  of  AS  is  1  m.  The  typical  range  cf  radar  target 
distances  is  a  few  to  hundreds  of  km.  The  required  aperture  size 
is 

L  =  XR/AS 

which  at  3  cm  wavelength,  a  target  distance  of  3  km  and  for  a 
required  resolution  of  1  m  evaluates  to  30  m.  In  this  example  the 
resolution  cell  is  30  times  smaller  than  the  aperture;  hence  the 
target  is  very  deep  in  the  near  field  (i.e.,  close  to  the 
antenna).  In  the  space-borne  illustration  of  Figure  1.4  the  wave- 
front  curvature  is  even  more  severe.  There  the  relation  of  aper¬ 
ture  size  to  cell  size  is  10  km  to  2  m  or  5000:1. 

The  practical  significance  of  near-field  operation  is  that  the 
imaging  system  must  focus  to  the  target  distance.  The  radio 
camera  automatically  performs  this  function  in  the  process  of 
adaptive  beamforming.  ABF,  being  designed  to  self-calibrate  the 
system  in  the  presence  of  arbitrary  phase  errors  across  the  array, 
easily  accommodates  for  the  quadratic  phase  variation  induced 
across  the  aperture  by  a  near-field  target.  The  two  groups  of  ABF 
techniques  that  have  been  developed  at  VFRC  compensate  for  aper¬ 
ture  distortion  and  also  focus  the  array  at  the  target  distance. 

A  distinction  can  be  made  between  self-calibrating  procedures 
according  to  whether  they  are  local  or  global.  A  local  procedure 
applies  adaptive  control  to  portions  within  the  entire  system  to 
correct  their  errors,  but  fails  to  compensate  for  all  the  errors 
in  the  system.  An  example  of  local  adaptivity  is  the  use  of  a 
pilot  signal,  i.e.,  a  known  waveform  at  a  known  time,  injected 
into  the  input  of  the  RF  section  of  a  phased  array.  By  measuring 
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the  output  signal  it  is  possible  to  deduce  the  distortions  in  the 
receiving  path  and  to  transmit  feedback  control  signals  to  the 
phase  shifters  and  amplifiers  for  the  purpose  of  compensating  for 
these  errors.  This  procedure  will  self-calibrate  the  receiver  but 
will  not  correct  for  element  position  errors  or  the  distortions 
caused  by  weather  or  propagation  conditions. 

Global  adaptivity  applies  feedback  control  over  the  entire  sys¬ 
tem,  simultaneously  correcting  for  propagation  medium  effects, 
geometric  distortion  in  the  aperture  and  phase  errors  in  the  dis¬ 
tribution  network. 

Adaptive  beamforming  as  applied  to  the  radio  camera  is  a  glob¬ 
ally  adaptive  procedure.  Instead  of  introducing  an  artificial 
pilot  signal  into  the  receiver,  as  is  the  practice  when  local 
adaptivity  is  to  be  achieved,  echoes  from  the  target  are  used  to 
derive  the  control  for  the  phase  shifters  or  the  weight  vector  of 
the  array.  These  echoes  experience  all  the  distortions  mentioned 
above.  They  interact  with  medium  turbulence  and  weather.  Their 
phases  are  distorted  by  element  position  error.  And  they  experi¬ 
ence  all  the  phase  and  amplitude  errors  due  to  circuit  mistunings, 
component  aging,  oscillator  drift,  etc.  It  is  these  distorted  sig¬ 
nals  that  are  converted  into  diffraction-limited  images  by  ABF.  To 
do  so  the  system  must  correct  all  errors  encountered  by  the 
echoes,  not  merely  the  ones  induced  in  the  receiver.  Consequently, 
the  ABF  control  procedure  is  global. 

A  companion  procedure  to  ABF  is  called  self-survey.  The  func¬ 
tion  of  a  self-survey  procedure  is  to  determine  the  locations  of 
the  parts  of  a  large  antenna  system,  for  once  the  element  posi- 
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tions  are  known,  it  is  then  possible  to  calculate  the  phase  shifts 
required  to  form  and  point  a  beam  in  an  arbitrary  direction.  A 
self-survey  system  employs  a  minimum  of  three,  and  typically  four 
to  six,  beacons  or  transmitters.  These  sources  are  distributed  in 
azimuth  throughout  the  sector  that  the  array  is  designed  to  scan. 
Either  time  of  arrival  or  phase  measurements  of  the  beacon  signals 
are  made  at  each  of  the  antenna  elements.  The  number  of  measure¬ 
ments  is  sufficient  to  calculate  the  entire  geometry  of  the  sys¬ 
tem,  i.e.,  the  positions  of  all  beacons  and  antenna  elements. 

Thus  a  two-step  operation  is  involved:  first  is  the  survey  of 
the  system  which  is  followed  by  beamforming  and  scanning.  ABF  is  a 
more  efficient  procedure.  It  recognizes  that  the  function  of  the 
antenna  array  is  only  the  latter  of  the  two  operations  performed 
by  the  self-survey  technique.  It  bypasses  the  survey  stage  com¬ 
pletely.  It  ignores  where  the  elements  are  located  and  therefore 
eliminates  the  requirement  to  calculate  their  locations.  It  only 
makes  those  measurements  that  permit  the  system  to  adjust  the 
phase  shifts  at  each  element  to  compensate  for  the  position  error. 
Furthermore,  the  requirement  of  the  self-survey  method  that  sever¬ 
al  beacons  be  deployed  over  the  scanning  sector,  which  limits  the 
use  of  the  procedure  to  fixed  installations  in  which  the  user  has 
control  of  the  nearby  terrain,  is  avoided  completely.  Instead,  in 
the  first  of  the  two  ABF  procedures,  only  a  single  source  in  the 
general  direction  in  which  the  system  is  scanning  is  need.  And 
experience  shows  that  this  source  does  not  have  to  be  a  beacon  or 
transmitter;  in  general,  some  scatterer  in  the  field  of  view  acts 
as  a  suitable  replacement.  Such  a  scatterer  is  called  a  beamform- 
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ing  source.  Even  that  modest  requirement  is  eliminated  in  the 
second  of  the  two  ABF  methods. 

There  are  several  desirable  properties  of  an  adaptive  beamform¬ 
ing  algorithm: 

1.  Independent  of  array  geometry  and  distortion.  The  algorithm 
should  be  unconstrained  by  the  locations  of  the  elements  or  the 
errors  in  their  positions.  In  simpler  terms,  the  self-calibration 
procedure  should  operate  successfully  without  regard  to  array 
design  or  the  nature  of  the  errors  in  the  system.  This  is  accom¬ 
plished  to  a  large  extent  in  the  algorithms  described  in  the 
following  sections. 

2.  Independent  of  source  distribution.  This  highly  desirable 
property  is  difficult  to  achieve.  The  algorithms  described  in 
Section  3  impose  conditions  upon  the  source  distribution  while 
those  in  Section  4  do  not. 

3.  Passive.  Adaptive  beamforming  is  most  desirable  when 
applied  in  the  passive  mode,  which  means  in  the  receiving  system 
only.  When  self-calibration  is  accomplished  in  this  manner,  the 
transmitter  and  receiver  functions  are  independent  and  the  equip¬ 
ment  need  not  be  colocated.  This  permits  bistatic  operation 
(receiver  and  transmitter  do  not  share  a  common  aperture)  and  pos¬ 
sibly  the  use  of  available  on-the-air  radiation  to  illuminate  the 
target  scene.  All  the  algorithms  have  this  desirable  property. 

4 .  All  algorithms  require  calculations  made  from  complex  field 
measurements  at  the  elements.  To  perform  a  calculation  based  upon 
all  the  measurements,  there  must  exist  within  the  system  an  intra¬ 
array  communication  subsystem  to  deliver  the  data  to  a  common 
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point  for  computing.  The  preferred  algorithms  are  those  that  mini¬ 
mize  this  additional  burden  upon  an  already  complex  system. 

5.  Noniterative.  The  calculations  referred  to  above  can  be 
made  at  the  element  level  or  aperture  plane  of  the  array,  which  is 
prior  to  the  image  processor,  or  in  the  output  or  image  plane  fol¬ 
lowing  the  image  processor.  Corrections  for  the  errors  are  gener¬ 
ally  applied  to  the  signals  entering  the  image  processor.  Thus  if 
the  calculations  are  made  in  the  aperture  plane,  they  can  be 
applied  directly  to  the  distorted  signals.  On  the  other  hand,  if 
they  are  derived  from  the  output  of  th'  image  processor  and  deliv¬ 
ered  to  the  input  signals  for  correction,  a  new  image  must  be 
formed  before  the  corrections  can  be  updated  or  a  new  weight  vec¬ 
tor  formed.  Such  an  iterative  procedure  across  the  image  processor 
is  costly  in  time  and  in  the  total  amount  of  computation  required. 
In  general,  a  noniterative  procedure  is  preferred. 

6.  Amplitude  only.  It  would  be  highly  desirable  if  amplitude 
measurements  would  suffice,  for  then  the  imaging  system  would  sim¬ 
plify  enormously.  Detection  would  be  by  an  envelope  detector  at 
each  element.  Without  phase  detectors  there  would  be  no  need  for 
frequency-  and  phase-locked  local  oscillators.  And  without  such 
LOs  there  would  be  no  need  for  an  LO  distribution  system  within 
the  array.  (See  [4],  chapters  5  and  8,  for  the  system  complexities 
that  result  from  the  need  to  preserve  phase.) 

No  procedure  has  been  discovered  as  yet  that  either  works  with¬ 
out  phase  information  or  that  reconstructs  phase  from  amplitude 
measurements  made  in  a  very  sparse  array.  However,  because  long 
distance  imaging  is  a  near-field  procedure  in  principle,  there  is 
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information  in  the  modulus  of  the  complex  field  measurement  that 
might  compensate  for  the  missing  phase  information.  The  reason  for 
his  conjecture  is  as  follows:  When  the  target  is  in  the  far  field 
the  relative  amplitudes  of  the  element  voltages  in  the  array  are 
independent  of  target  distance.  Not  so  in  the  near  field,  however, 
for  there  the  phase  and  amplitude  of  the  radiation  field  couple. 
This  means  that  the  set  of  amplitudes  across  the  array  due  to 
radiation  or  reflection  from  a  target  scene  that  is  located  in  the 
near  field  is  a  function  of  distance  to  the  target.  Consequently, 
amplitude  measurements  made  at  two  or  more  distances  from  the  tar¬ 
get  convey  information  about  the  scene  that  is  normally  carried 
only  in  the  arguments  of  the  measurements. 

The  second  possibility  -  the  reconstruction  of  phase  from 
amplitude  measurements  -  is  a  subject  of  current  research.  As  this 
field  matures,  solutions  may  arise  useful  to  the  large  array  prob¬ 
lem.  One  promising  technique  is  reported  by  Tsao.  Only  10%  of  the 
receivers  in  the  array  require  phase  measurement;  90%  of  the  array 
consists  of  amplitude-only  detectors  [7]. 
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3.  DOMINANT  SCATTERER  ALGORITHMS  (DSA) 

3.1  ADAPTIVE  BEAMFORMING  ON  A  SINGLE  DOMINANT  SCATTERER 

The  primary  procedure  used  at  VFRC  today  assumes  that  there 
exists  somewhere  in  the  field  of  view  (FOV)  of  the  imaging  system 
a  point-like  scatterer  or  source  having  large  radar  cross-section 
or  source  strength.  The  theory  governing  the  requirement  on  cross- 
section  or  source  strength  is  given  in  [8]  .  Such  a  target  or 

source  radiates  a  nearly  spherical  wavefront,  thereby  inducing  a 

simple  phase  variation  across  the  array.  Provided  that  the  array 
is  linear  or  planar,  the  phase  variation  is  linear  if  the  target 
or  the  source  is  in  the  far  field,  or  approximately  quadratic  when 
the  source  is  in  the  near  field.  Deviations  from  such  simple 

behavior  indicate  geometric  distortion  in  the  array,  electrical 

mistunings  or  wavefront  distortion  due  to  turbulence  in  the  propa¬ 
gation  medium. 

Consider  a  one-dimensional,  far-field  source  distribution  s  = 
s(u) .  Its  radiation  field  in  the  axis  of  the  array  is 

e  =  e  (x)  =  F-1  {s  (u)  } 

where  the  right  hand  side  is  the  inverse  Fourier  transform  of  the 
source  distribution  s(u),  and  u  =  sin0  where  0  is  the  angle  from 

the  normal  to  the  array.  For  simplicity  operation  notation  is  used 
in  the  following  development.  F  indicates  the  Fourier  transform 
and  F-1  its  inverse.  The  electric  field  in  the  x-axis  is  then 
written 

e  =  e (x)  =  F_1s  (3.1) 
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Now  imagine  that  the  system  could  measure  e  and  form  a  weight  vec¬ 
tor  w  equal  to  its  conjugate: 

w  =  w(x)  =  e*  =  Fs*  (3.2) 

The  radiation  pattern  of  the  system  is  the  inverse  transform  of  w, 
which  is 

f  =  f ( u )  =  F_1w  =  F-1Fs*  =  s*  (3.3) 

Thus  the  radiation  pattern  is  the  complex  conjugate  of  the  source 
distribution1.  The  most  interesting  aspect  of  this  procedure  is 
that  no  account  was  taken  of  the  shape  of  the  antenna,  the  place¬ 
ment  of  its  parts  or  geometric  or  electrical  errors  in  the 
antenna.  To  utilize  this  relationship  in  an  ABF  algorithm,  whereby 
the  shape  of  the  radiation  pattern  is  the  modulus  of  the  angular 
profile  of  the  source,  it  is  necessary  that  there  be  at  least  one 
scatterer  of  small  physical  size  that  has  a  radar  cross  section 
(RCS)  so  large  that  its  reflection  dominates  all  the  other  scat- 
terers  in  its  range  cell2.  A  corner  reflector  is  an  example;  its 
RCS  is  approximately  T4/X2  where  T  is  the  length  of  one  side  and 
its  reflection  profile  approximates  a  5-function.  Let  this  reflec¬ 
tor  be  at  distance  r  =  rQ .  The  source  distribution  at  this  range 
is 

sQ  =  s0(u)  »  5(u-u0)  (3.4) 

where  u0  is  the  bearing  to  the  dominant  scatterer. 

xThe  weight  vector  is  chosen  as  the  complex  conjugate  of  the  electric  field  to 
ensure  that  the  radiation  pattern,  when  formed  from  the  e-field  of  a  point 
source,  points  to  that  source.  This  is  called  ret rodi recti ve  beamforming.  If 
the  conjugate  is  not  employed  the  beam  points  in  the  mirror  or  Snells-law 
direction. 

2A  range  cell  or  range  bin  is  an  interval  in  range  such  that  all  scatterers 
within  it  return  echoes  that  overlap  in  time  in  the  receiver.  Its  nominal 
extent  is  ct/2  where  c  is  the  speed  of  light  and  X  is  the  radar  pulse 
duration. 
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Next,  let  the  adaptive  system  recognize  the  presence  of  the 
single,  dominant  scatterer  at  r  =  rQ  (by  the  procedure  described 

in  Section  3.2)  and  form  its  weight  vector  from  measurements  of 
the  radiation  from  that  range.  Then  from  (3.3)  and  (3.4) 

f  =  f (u)  =  8(u-u0)  (3.5) 

is  the  radiation  pattern  of  the  antenna.  When  this  radiation  pat¬ 
tern  scans  the  source  distribution  s^  =  s^(u)  at  range  Ri,  the 

image  delivered  by  the  system  is 

s =  s^(u)  =  Jf  (x)  si  (u-x)  dx  £  f(u)*s^(u) 

«  5  (U-UQ)  *s.^  (u) 

*  Si(u-u0)  (3.6) 

Thus  the  output  is  an  excellent  replica  of  the  target  scene  except 
for  the  unknown  displacement  angular  uQ. 

We  now  introduce  an  antenna  into  the  equations.  Let  the 

antenna  be  a  phased  array  which  is  designed  as  a  sampled  aperture 

of  finite  length.  Let  the  design  location  of  the  nth  element  be  xon 

and  the  weight  won.  The  aperture,  as  designed,  can  be  described  by 

N 

wo  =  wo(x>  “  X  won5(x_xon>  <3-7> 

where  the  weights,  as  a  rule,  are  real,  and  are  selected  to  pro¬ 
vide  shaping  to  the  side  radiation  pattern.  Because  of  gain  errors 
and  aperture  distortion,  the  actual  weight  vector  differs  from 
(3.7)  and  can  be  written 

w(x)  =  I  wn8(x-xn)  (3.8) 

where 


w 


n 


*=  ane^n 


(3.9) 
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Now  the  weights  have  become  complex.  The  wn  differ  from  the  design 
values  won  primarily  due  to  the  phase  errors  in  the  system  and  the 
xn  differ  from  the  xon  due  to  the  geometric  errors  in  the  locations 
of  the  elements  of  the  antenna.  For  the  most  part  the  modulii  of 
the  wn  do  not  differ  significantly  from  the  design  values  because 
the  wavefront  of  the  radiation  from  a  point  source  is  spherical; 
hence  a^  **  ^on  * 

The  weight  vector  can  be  viewed  as  a  transducer  that  converts 
the  electric  field  across  the  aperture  to  a  current  i  =  i(x)  that 
drives  the  electronic  system.  i  is  given  by 

i  =  we  =  wF-1s  (3.10) 

Let  the  adaptive  beamforming  subsystem  measure  i  and  form  a  com¬ 
pensation  or  correction  weight  vector 

c  =  c(x)  =  i*  *  w*e*  =  w*Fs*  (3.11) 

The  composite  weight  vector  is  the  cascade  or  product 

wc  =  ww*Fs*  (3.12) 

The  term 

WW  =  I  X  wnwm  5(x-xn)8(x-xm) 
n  m 

*  2  A 

=  X  wnwn  8(x-xn)  =  X  an  5(x-xn)  _  wx  (3.13) 

n  n 

is  a  sampling  function  having  sample  points  at  the  actual  element 

2 

location  xn,  and  real  weights  an,  implying  that  the  phase  errors 

in  the  aperture  have  been  eliminated.  The  composite  weight  vector 
wc  =  w^s*  (3.14) 

and  the  radiation  pattern  is 
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f  =  f(u)  =  F“^wc  =  F  ^v^Fs*  =  F“1w1  *  s*  (3.15) 

where  F”^w1  =  (u)  is  the  radiation  pattern  of  the  sampled  aper¬ 

ture  defined  by  (3.13).  Assuming,  as  before,  a  source  distribution 
at  r  =  rQ  given  by  (3.4)  and  the  composite  weight  vector  (3.12) 

formed  from  measurements  of  the  echoes  from  that  range,  the  adap¬ 
tively  formed  radiation  pattern  is 

f(u)  *  fj (u)  *  8(u-u0)  =  f1(u-uQ)  (3.16) 

The  key  feature  of  (3.16)  is  that  all  phase  errors  have  been 

eliminated  by  the  adaptively-formed  compensation  weight-vector. 

2 

The  aperture  weights  an  are  real,  as  they  should  be.  They  differ 

from  the  design  values  won .  However,  because  of  the  approximate 
relationship  an  =*  won,  that  difference  can  be  easily  corrected  by 

modifying  the  initial  won  to  V won  to  account  for  the  fact  that  the 

2 

final  weights  are  an.  Alternately  the  correction  can  be  introduced 

in  the  signal  processing  computer.  The  simplest  procedure  (used 
for  all  VFRC  imaging)  is  to  set  won  =  1  for  all  n,  and  to  perform 
the  desired  aperture  tapering  upon  the  stored  data  records  in  the 
computer.  With  this  procedure  the  square-root  relationship  between 
the  weights  of  wx  and  wQ  is  avoided  and  an  arbitrary  taper  most 

suitable  to  the  particular  imaging  problem  can  be  applied. 

3.2  THE  RADIO  CAMERA 

The  adaptive,  retrodirect ive  beamforming  procedure  described 
above  is  the  key  to  the  design  of  the  radio  camera.  The  radio  cam¬ 
era  is  a  large  imaging  system  that  operates  at  wavelengths  very 
large  compared  to  optics,  and  its  large  size  implies  that  it 
suffers  from  distortions  in  its  shape  and  in  errors  in  the  loca- 
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tions  of  its  parts.  Propagation  and  circuit  errors  are  assumed  as 
well.  The  means  for  applying  the  procedure  given  in  the  preceding 
section  is  described  in  detail  in  Chapter  5  of  [4]  and  in  [3]. 
Figure  3.1  illustrates  the  logic  by  which  calibration  information 
is  extracted  from  the  radiation  field  from  a  target.  A  point 
source  or  target,  identified  as  a  phase  synchronizing  source  in 
the  sketch,  reflects  a  transmitted  pulse  toward  the  array.  Because 
of  the  small  size  of  the  target,  it  reradiates  a  spherical  radia¬ 
tion  field,  meaning  that  the  amplitude  and  phase  are  constant  on 
any  sphere  centered  on  the  synchronizing  source.  Assuming  free- 
space  radiation,  the  complex  field  amplitudes  measured  by  the  ele¬ 
ments  of,  say,  a  linear  or  planar  array,  would  have  equal  moduli. 
Their  phases,  however,  would  be  random  and  meaningless  because  of 
the  geometrical  distortion,  even  with  element  position  errors  as 
little  as  several  tenths  of  a  wavelength.  Any  circuit  mistunings 
or  mismatches  would  further  distort  the  measured  phases.  The  sys¬ 
tem  compensates  for  the  phase  distortion  across  the  aperture  by 
comparing  the  phases  measured  at  the  various  elements  and  deriving 
feedback-control  signals  from  them  to  drive  the  first  bank  of 
phase  shifters  so  that  the  signals  passing  through  the  phase 
shifters  are  all  cophased  (described  below) .  Doing  so  phase  shifts 
each  receiver  channel  by  the  negative  of  the  received  signal. 
Because  the  amplitudes  are  equal  and  the  phase  shifts  are  the  neg¬ 
atives  of  the  measured  signal  phases,  the  weight  vector  is  the 
complex  conjugate  of  the  field  measurements.  The  first  bank  of 
phase  shifters  and  the  circuit  that  controls  it  is  the  calibration 
or  compensation  subsystem.  The  sum  of  the  cophased  signals  is 
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FIGURE  s,  1.  Distorted  array  self -coheres  on  a  phase- synchronizing  source:  >  -  antenna  ele¬ 
ments;  •>  -  reference  element;  L  -  size  of  array;  R  -  nominal  distance  to  array;  r„  -  distance  to 
nth  element  (From  [4].) 
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that  of  an  array  focused  on  the  phase  synchronizing  target.  The 

two-dimensional  radiation  pattern  of  the  array  is  centered  at  this 

Xr 

point.  The  beam  cross  section  in  the  focal  zone  is  ~  in  the  lat- 

Li 

eral  or  cross-range  dimension  and  approximately  7X(R/L)2  in  the 
range  or  longitudinal  dimension. 

Once  the  beam  is  formed,  a  second  set  of  phase  shifters  can 
scan  it  to  form  a  map  or  an  image  of  a  target.  In  essence,  the  in¬ 
puts  to  the  second  bank  of  phase  shifters  can  be  viewed  as  the 
outputs  of  antenna  elements  of  a  properly  designed,  error-free 
array.  Scanning  is  an  open-loop  procedure,  exactly  as  in  a  conven¬ 
tional  phased  array.  The  extent  to  which  scanning  can  be  accom¬ 
plished  is  limited,  however,  by  the  amount  of  geometrical  distor¬ 
tion  in  the  array.  This  is  because  only  a  single  beamforming 
source  is  used  to  calibrate  the  system.  The  calibration  is  correct 
in  the  direction  of  the  beamforming  source  and  at  its  distance. 
Residual  distortions  develop  as  the  beam  is  scanned  away  from  the 
location  of  the  beamformer.  Errors  in  the  array  in  the  direction 
of  the  beamforming  source  and  at  its  range  are  fully  compensated. 
The  errors  which  dominate  under  scanning  are  the  distortions  in 
the  axis  of  the  array  [4,  Section  7.3].  (Experiments  supporting 
the  theory  are  found  in  the  same  reference.)  Errors  normal  to  the 
array  direction  and  to  the  direction  to  the  beamforming  source 
have  a  relatively  minor  effect. 

The  tolerance  theory  for  scanning  in  the  near  field  of  the 
array  is  given  in  [4,  Section  7.5].  There  it  is  shown  that  the 
maximum  scan  angle  from  the  focal  point  is 
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e 


max 


4KGy 


(3 . 17; 


where  Cx  is  the  standard  deviation  of  a  zero-mean  random  variable 

in  the  x-position  of  each  of  the  antenna  elements  in  the  array. 
Equation  (3.17)  is  derived  on  the  assumption  that  the  maximum  nom¬ 
inal  tolerance  on  the  loss  of  array  gain  is  1  dB .  The  angular 
field  of  view  is  twice  0max  given  by  (3.17)  and  the  cross-range 

dimension  of  the  sector  that  may  be  scanned  by  the  adaptively 


formed  beam  is  R  times  this  value,  as  given  in  (3.18) 

Xr 

AS  -  = 


'max1 


2  no 


(3.18) 


X 


The  extent  in  range  AP  over  which  the  beam  may  be  scanned  is  much 
larger.  It  is  shown  in  [4,  Section  7.5]  that  the  ratio  of  the 


dimensions  of  the  scanning  sector  is  given  by 
45.  =  Vl2  R 
AS  L 


(3.19) 


where  R  is  the  range  to  the  focal  zone  and  L  is  the  size  of  the 
array.  In  general,  AR  »  AS  because  normally  R  »  L.  Thus  the 

position-error  limitation  to  scanning  is  most  heavily  felt  in  the 
cross-range  or  angle  direction. 

The  two  banks  of  phase  shifters  in  Figure  3.1  indicate  only 
that  two  distinct  functions  must  be  performed.  These  functions  can 
be  performed  in  hardware  or  within  a  computer.  The  procedure  used 
at  VFRC  is  illustrated  in  Figure  3.2.  A  search  procedure  is  indi¬ 
cated  by  which  the  imaging  system  discovers  the  phase  synchroniz¬ 
ing  source,  which  is  a  critical  operation  for  the  success  of  the 
system. 
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FIGURE3-2.  Radio  camera  data  format  and  procedures.  (From  B.  D.  Steinberg,  IEEE  Trans. 
Antennas  and  Propagation,  Ap-29,  No.  5,  September  1981.) 
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The  search  procedure  is  based  upon  a  minimum  amplitude  vari¬ 
ance  (MAV)  test  upon  the  echoes  streaming  across  the  array. 
Referring  to  Figure  3.2,  a  transmitter  launches  a  pulse  toward  a 
target  area,  in  which  is  embedded  a  small  target  of  large  radar 
cross  section.  A  corner  reflector  is  illustrated  in  the  figure, 
although  experience  shows  that  it  is  rarely  needed  because  targets 
of  opportunity  having  appropriate  scattering  properties  are  plen¬ 
tiful  (Section  3.5).  The  reflected  echo  from  the  range  bin  con¬ 
taining  the  reference  target  is  dominated  by  the  reflection  from 
that  target.  The  distorted  array  consists  of  antenna  elements 
followed  by  receivers  that  coherently  detect  and  demodulate  the 
signals.  The  data  are  sampled  and  stored  as  shown  in  the  lower 
right  of  the  figure.  The  sampling  interval  is  approximately  the 
pulse  duration.  Each  row  of  data  is  the  echo  trace  received  by  a 
single  element.  Successive  rows  correspond  to  successive  elements 
in  the  array.  Each  column  of  data  is  called  a  range  bin.  The  data 
represent  a  snapshot  of  the  radiation  field  from  sources  at  a  dis¬ 
tance  corresponding  to  the  time  of  arrival  of  those  signals.  For 
simplicity,  it  is  assumed  that  range  walk,  or  range-bin  slippage 
has  been  corrected. 

To  locate  the  range  bin  in  which  the  target  most  closely 
resembles  a  point  source,  the  signal  processor  calculates  the  nor¬ 
malized  variance  of  the  measured  amplitudes  of  the  signals  re¬ 
ceived  across  the  array  in  each  range  bin.  The  range  bin  with  the 
smallest  variance  is  chosen,  for  that  echo  sequence  most  closely 

approximates  the  expected  field  from  a  point  source.  Designate 
that  range  by  R0 .  The  weight  vector  adopted  by  the  system  is  the 
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complex  conjugate  of  the  complex  column  vector  of  signals  received 
at  that  range.  It  is  shown  in  Chapter  6  of  [4],  and  also  in  [8], 
that  the  reference  target  echo  need  be  no  more  than  4  dB  greater 
than  noise  and  clutter  in  that  range  bin  for  successful  phase  syn¬ 
chronization  of  the  array  to  be  achieved. 

Table  3.1  describes  the  algorithm.  The  measured  signals  from 
the  ith  range  bin  at  the  nth  element,  corrected  by  the  shape  of 
the  element  pattern,  are  Aj.nexp  ( j\yin)  .  Step  3  finds  the  reference 
range  RQ  in  which  the  amplitudes  Aon  are  nearly  equal.  Because  the 

Aon  are  nearly  equal,  the  measured  value  at  the  nth  element  can  be 
written  approximately  as  A  exp(j\|Ton).  Its  conjugate  A  exp(-j\|/0n)  is 

the  weight  vector.  Step  4  gives  the  phase-corrected  signals  in  the 
Oth  range  bin.  The  corrected  inth  term  is  proportional  to 
Ainexp[  j  (Von-Von)  ]  (Step  5)  . 

This  term  has  now  been  corrected  for  array  distortion.  It 

also  has  been  corrected  for  the  near-field  curvature  of  the  phase- 
front  associated  with  the  distance  RQ  to  the  phase  synchronizing 

source.  Step  6  alters  the  phase  curvature  by  replacing  the 
quadratic  phase  variation  of  the  reference  range  by  the  proper 
curvature  for  the  ith  range  bin,  and  similarly  for  all  other  bins. 
Following  this  step,  the  beam  may  be  scanned  in  angle,  for  it  is 
now  in  focus  at  all  ranges.  Steps  7  and  8  complete  this  process. 

High  resolution  imaging  implies  that  the  resolution  cell  in 
the  transverse  dimension  is  smaller  than  the  array.  The  target 
always  is  in  the  near  field  and  the  array  must  therefore  be 
focused.  Adaptive  beamforming  automatically  focuses  the  array  to 
RQ  (step  5)  .  However,  it  is  rarely  necessary  to  perform  step  6 
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Step 

1 

2 

3 

4 

5 

6 
7 


TABLE  3-L  Steps  in  Radio  Camera  Imaging  [ 4  ] 


Measure  and  store  complex  en¬ 
velopes  of  echo  samples 

Correct  amplitudes  by  dividing 
by  element  pattern  estimate  /„ 
Find  R0  such  that  A0„  •  A, 
all  n 

Phase  rotate  at  R0  by  phase 
conjugate  in  relation  to  refer¬ 
ence  element,  exp/(tfa>  -  i£0n) 
Phase  rotate  at  all  range  ele¬ 
ments 

Focus  at  each  range  R, 

Phase  shift  linearly  with  angle 


V  e1*- 

w  tar 

range  bin—'*  ^-element  number 
A  e*'* 

Ae 


A  iKeJ(*‘>~  *«•**<») 

A  l/«o)]  A  R 

B,.e~ikx-“ 


*,(«)-  E  BiHe-Jkx- 

It—  I 


8 


Sum  at  each  range  element 


which  refocuses  the  array  simultaneously  to  all  range  bins.  At 
most  it  is  necessary  to  refocus  in  steps  of  100s  of  meters.  This 
is  because  the  depth  of  field  of  an  array  usually  is  larger  than 
the  extent  in  the  range  dimension  of  the  sector  being  imaged. 

The  depth  of  field  of  an  aperture  (array,  lens,  dish,  etc.) 
is  the  extent  in  range  over  which  the  aperture  is  nominally  in 
focus  and  is  approximately  7A.(R/L)2  [9].  The  cross-range  resolution 

is  approximately  AS  =  Xr/L.  Combining, 

R  o  AS  o  7  o 

d.o.f.  =  iX  (-)  =  iX  (—)  -  r  (As)2. 

Jj  A,  A. 

Given  a  3-cm  system  focused  to  a  1-m  beam  cross  section,  the  depth 

of  field  is  about  200  meters. 

Figure  3.3  (top)  illustrates  the  spatial  distribution  of 
scatterers  of  various  sizes  located  at  some  range  R.  The  illus¬ 
tration  to  the  left  shows  a  group  of  small  scatterers  of  random 
sizes  and  locations.  This  distribution  models  statistically  homo¬ 
geneous  clutter,  uniformly  distributed  throughout  the  illuminated 
sector.  On  the  right  one  large  target  or  scatterer  has  been 
added.  The  bottom  sketches  picture  the  difference  between  the 
phasefronts  of  the  fields  received  at  the  array  for  the  two  cases. 
The  phasefront  on  the  left  is  a  sample  of  a  zero-mean  random  pro¬ 
cess  whereas  the  phasefront  on  the  right  has  a  well  defined  slope 
determined  by  the  bearing  from  the  array  to  the  large  scatterer. 
Superimposed  on  the  average  phasefront  is  a  perturbation  field 
caused  by  the  radiation  from  all  the  remaining,  small  scatterers 
in  the  range  bin. 
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Figure  3 . 3  Phase  distributions  across  array.  From  SteinbeTg,  J.  Franklin  Inst.  (December 
1973). 
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The  argument  that  follows  is  based  upon  plane  wave  propaga¬ 
tion,  but  it  is  equally  valid  in  the  near  field.  The  sum  of  the 
waves  emanating  from  the  scatterers  is  the  radiation  field  at  the 
receiving  array.  Only  those  scatterers  illuminated  by  the  trans¬ 
mitting  beam  contribute  to  the  reradiation  field.  Call  the  source 
function  of  this  field  p(u) .  We  can  write  this  source  function 
p  =  p(u)  =  fT(u)  s(u)  (3.20) 

fT  =  fT(u)  is  the  radiation  pattern  of  the  transmitter  and 

s  =  s  (u)  =  X  SjL8(u-Uj_)  (3.21) 

i 

is  the  scatterer  distribution,  Sj_  and  u^  being  the  complex 

strength  and  bearing  of  the  ith  scatterer.  The  electric  field  in 

the  receiving  array  is 

e  *  e  (x)  *  F-1p 

=  F-1fT*F-1  X  si5(u-ui) 
i 

=  wT*  X  s^exp ( jkxu^)  (3.22) 

i 

where  wT  is  the  aperture  excitation  of  the  transmitter  and  F_1 

means  inverse  Fourier  transform.  The  electric  field  is  the  sum  of 
randomly  weighted  (by  the  s^)  plane  waves  arriving  from  random 

directions  u^,  each  convolved  with  wT. 

Let  1  be  the  size  of  the  transmitting  antenna  and  L  the  size 
of  the  receiving  antenna.  The  transmitted  beamwidth  A0T  *  X/l  and 

the  width  of  the  clutter  patch  illuminated  at  range  R  by  the  beam 
is  As  =  XR/1.  The  reradiation  is  lobular  with  typical  lobe  width 
A,/As  =  1/R.  At  the  receiver,  also  a  distance  R  from  the  clutter 
patch,  the  lobe  cross  section  is  (1/R)R  =  1,  which  is  approxi- 
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mately  the  separation  between  such  lobes.  These  L/l  lobes  across 
the  receiving  array  are  statistically  independent  because  they  are 
derived  from  the  Fourier  transform  of  a  band-limited  random  pro¬ 
cess.  Thus  the  phase  of  the  e-field  across  the  receiving  array  is 
a  band-limited  random  variable  characterized  by  the  order  of  L/l 
random  values  in  the  interval  [-71,  Ji] .  Such  a  field  is  useless  for 

phase  synchronizing  the  receiving  array. 

The  situation  changes  dramatically  when  a  large  radar  cross 
section  (RCS)  reflector  is  added  (upper  right) .  The  sum  of  the 
scattered  fields  from  the  clutter  does  not  change,  but  added  to  it 
is  the  strong  plane  wave  from  the  prominent  scatterer.  The  e-field 
can  be  written 

e  =  e(x)  =  s0exp(jkxu0)  +  X  s^exp ( jkxu^)  (3.23) 

i 

in  which  the  first  term  dominates.  This  field  can  be  viewed  as  a 
perturbed  plane  wave,  i.e.,  its  direction  of  arrival  is  approxi¬ 
mately  uQ,  meaning  that  its  phase  slope  is  nearly  but  not  exactly 

kuQ  at  every  point.  The  phase  of  the  e-field  (3.23)  is 

s0sinkxu0  +  2  si  sinkxui 

<J)X  =  tan*1  - 4 - 

s0coskxu0  +  X  si  coskxui 
i 

X  Si  sin  kx(ui-u0) 

=  kxuc  +  tan-1  [ - - -  ]  (3.24) 

s0  +  I  Si  coskx(ui~u0) 
i 

The  sum  field  from  the  random  scatterers  can  be  called 

A(x)expja(x)  =  X  siexp(jkxui)  (3.25) 

i 
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With  A  «  sQ  throughout  the  array  when  sQ  is  large  compared  to  the 

sum  of  all  the  other  scatterers, 

A(x)  sin  a(x) 

4>x  =  kxuQ  +  tan-1  {  - - -  } 

So 

A(x)  .  .  . 

»  kxuQ  +  — —  sm  a(x) 
u  s0 

=  kxuQ  +  &}>(x)  (3.26) 

where  5<|>(x)  is  the  deviation  of  <)>(x)  from  kxuQ. 

Equation  (3.26)  is  the  phase  available  to  be  measured  by  the 
system  for  self-calibration  of  the  array.  The  first  term  is  the 

linear  phase  term  with  constant  slope,  which  is  the  desired  term. 
The  second  term  is  the  error  phase.  The  dominant  field  due  to  sQ 

need  exceed  the  scattered  field  (3.25)  by  no  more  than  6  dB  to 
enable  the  dominant  scatterer  algorithm  (DSA)  to  extract  the  phase 
information  it  needs  with  adequate  accuracy  to  form  and  scan  a 
main  beam  [4,  Section  7.6]. 

The  procedure  is  summarized  in  the  mathematics  that  follows. 
Consider  a  linear  array  with  N  elements  located  at  (^n'^n^'  1  -  n 
<  N.  The  yn  are  the  deviations  in  the  element  positions  in  the 

array  that  are  to  be  corrected  by  the  adaptive  beamforming  proce¬ 
dure.  Let  P  ,  1  £  n  <i  N,  0  ^  IP  I  :£  K,  denote  the  electronic  phase 
n  n 

error  at  each  of  these  elements.  These  errors  are  caused  by  mis- 
tunings,  impedance  mismatches,  oscillator  drifts,  etc.  Assume  that 

somewhere  in  the  FOV  of  the  imaging  system  there  exists  a  range 
bin  containing  a  single  dominant  scatterer  plus  clutter.  Let  RQ 

denote  the  distance  of  this  range  bin  from  the  origin,  which  is 
assumed  to  be  located  at  the  first  element  (x^y^O)  (see  Figure 
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3.9).  Let  the  bearing  to  the  scatterer  be  uQ  =  sin0Q,  where  0Q  is 

the  angle  measured  clockwise  from  broadside  to  the  array. 

The  phase  of  the  echo  from  the  RQ  range  bin  measured  at  the 


nth  receiver  is 


0, 


o,n 


krn  n  +  B  +  &|>_ 
°'n  °,n 


(3.27) 


where  k  =  2it/Xf  rQ^n  is  the  distance  from  the  beamforming  scat¬ 
terer  (also  called  adaptive  beamforming  source,  beamformer  or  ref¬ 
erence  target)  to  the  nth  receiver,  and  is  equal  to 


ro,n  t (xn  ~  Rosin0o) 2  +  (yn  -  Rocos0o)^] 


1/2 


(3.28) 


Pn  is  the  electronic  phase  error  at  the  nth  receiver  and  6<))0  n  is 


the  phase  perturbation  in  the  spherical  wavefront  of  the  beamform¬ 
ing  scatterer  due  to  clutter  and  noise. 

The  Fresnel  approximation  for  the  distance  rQ,n  can  be  used 

in  most  microwave  imaging  scenarios.  Further,  the  perturbation  of 
the  nth  element  in  a  direction  perpendicular  to  the  array,  yn, 
will  be  much  smaller  than  its  location  along  the  array,  i.e.,  xn . 

Expanding  (3.28)  in  a  Taylor's  series  and  ignoring  all  terms  other 
than  the  linear  and  quadratic  terms  in  X  and  Y  yields 


ro,n  Ro  ”  xnuo  ~  ¥n  ^  uo^  ^  +  xn  ^  ”  un^Ro-  (3.2  9) 

A  similar  expression  pertains  to  any  other  range  bin.  Thus,  rm  n, 

the  distance  from  a  scatterer  in  the  mth  range  bin  to  the  nth  ele¬ 
ment,  is  found  from  (3.29)  with  all  'o'  subscripts  replaced  by 
'm'  . 
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The  signal  phase  <|>m  ^  at  the  array  origin  is  treated  as  a 


reference  phase  for  the  mth  bin  and  is  subtracted  from  all  the 

phases  measured  in  the  array,  thereby  removing  the  target  distance 
Rm.  In  the  reference  range  bin  RQ/  the  difference  phase  at  the  nth 


element  is  (substituting  (3.29)  into  (3.27)) 


2,1/2 


4>o,n  =  4>o,n  ~  4*0,1  =  k{xnuo  +  Vn*1  ~  uo} 

~  kxn(1  "  uo)  7  2Ro  +  Pn  "  bl  +  ^o,n 


(3.30) 


where  The  negative  (or  conjugate  phase)  of 

/  **  o,  n 


(3.30)  is  the  phase  shift  introduced  into  the  nth  channel  to  com¬ 
pensate  for  array  distortion  and  electrical  phase  errors  in  the 
receivers.  This  is  the  ABF  process. 

It  is  desirable  to  combine  the  second,  fourth  and  fifth  terms 
to  form  a  phase  error  (3n  =  |3n  -  (3-^  +  kyn(l-u0)  '  and  to  treat  it 


as  being  nearly  independent  of  the  viewing  angle  of  the  imaging 
system.  The  phase  errors  Pn  and  (3  1  are  electronic  errors  and 


therefore  are  independent  of  viewing  angle.  The  last  term,  on  the 
other  hand,  is  not;  it  is  an  explicit  function  of  the  angle,  which 
for  any  scanning  problem  can  vary  only  over  the  width  of  the  imag¬ 
ing  sector.  However, the  FOV  of  most  microwave  imaging  systems  is 
narrow  because  the  transmitter  beamwidth  is  typically  the  order  of 
a  few  degrees.  Hence  this  error  may  be  assumed  to  be  independent 

of  this  small  variation.  Consequently,  summing  the  electronic 
phase  error  (3  -  (3^  and  the  perturbation  error  kyn(l  -  u  )  '  gives 
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a  phase  error  Pn  that  is  almost  independent  of  the  beamformer 
direction  uQ.  Using  pn  =  (3n  -  +  kyn(l  ~  u^)1^2  in  (3.30)  yields 

2  2 
kxnuc  -  kxn(l  -  uQ) 

=  2Rq  +  ^n  +  ^o,n*  (3.31) 

The  first  two  terms  point  and  focus  the  array  to  (uQ,R0) .  The 
third  term,  pn,  is  the  phase  error  from  the  RQ  bin  caused  by  elec¬ 
trical  phase  distortion  and  element  position  errors.  Except  for 
the  small  variation  across  the  FOV,  this  error  is  constant  for  all 
range  bins;  hence  its  inclusion  in  the  weight  vector  of  the  array 
eliminates  the  primary  source  of  distortion.  The  fourth  term  is 
caused  by  clutter  and  noise.  Because  the  clutter-induced  phase 
perturbation  is  not  constant  from  range  bin  to  bin,  the  6<t>0^  n  term 

does  not  compensate  for  this  source  of  phase  error.  Let  5<J)m>n  be 

the  clutter-induced  phase  perturbation  in  the  mth  range  bin;  then 
5V  n  -  **o,n  is  the  residual  phase  error  after  adaptive  beamform¬ 
ing  and  is  responsible  for  the  loss  in  mainbeam  gain  and  the  rise 
in  the  sidelobe  level  that  result  from  the  ABF  process. 

The  nth  component  of  the  phase-compensating  weight  vector  is 
wn, DSA  "  e*PM<l>o,n>  <3-32) 

Thus  the  DSA  phase  corrections  compensate  for  both  electronic 
phase  errors  and  perturbation  of  the  receivers  along  the  Y-axis; 
however,  they  are  contaminated  by  clutter  and  focus  the  array  to 
uQ  and  not  to  the  center  of  the  FOV,  which  we  may  call  u.  Focusing 
the  array  to  uQ  leads  to  a  pointing  error  u-uQ  in  the  target  loca¬ 
tion  estimates.  The  multiple  scatterer  algorithm  (MSA),  described 
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in  Section  3.5,  reduces  the  clutter  noise  in  the  phase  corrections 
and  yields  a  better  estimate  of  the  target  location. 

3.3  DSA  EXPERIMENTS 

Radio  camera  images  have  been  made  with  the  DSA  applied  to 
data  from  three  experimental  facilities.  The  first  is  a  cable-sup¬ 
ported  synthetic  receiving  array  83  m  in  length.  This  array  is 
2500  wavelengths  at  the  3  cm  wavelength  at  which  it  was  operating. 
Its  resolving  power  is  comparable  to  human  vision.  The  experimen¬ 
tal  procedure  is  shown  schematically  in  Figure  3.4.  A  single  re¬ 
ceiver  was  moved  along  a  cable  10  m  above  the  ground.  The  cable 
was  83  m  long  and  supported  only  at  the  ends  and  center.  The  radar 
transmitter  was  directed  approximately  normal  to  the  array.  The 
receiver  with  1  ft  dish  was  pulled  from  end  to  end  in  approxi¬ 
mately  two  minutes.  It  was  turned  on  330  times  at  predetermined 
locations  as  it  moved  over  the  83  m  course.  The  estimated  uncer¬ 
tainty  in  antenna  location  at  each  of  the  sample  points  was  the 
order  of  1  wavelength  rms. 

Figure  3 . 4A  shows  images  made  with  this  array.  The  scene  is 
of  rows  of  houses  on  streets  in  the  town  of  Phoenixville  with  and 
without  adaptive  beamforming.  The  same  microwave  data  were  used  to 
form  these  images.  The  dominant  scatterer  adaptive  beamforming 
algorithm  was  used  to  form  the  upper  image.  Careful  examination  of 
the  image  data  shows  that  the  "corrected"  array  was  operating 
diffraction  limited.  The  lower  picture,  without  adaptive  beamform¬ 
ing,  is  a  useless  array  of  pixels.  The  image  of  the  housing  devel¬ 
opment  in  Figure  3.4B  was  also  made  with  this  array. 
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FIGURE  3.4  CABLE-SUPPORTED  SYNTHETIC  RECEIVING  APERTURE 
AT  THE  VALLEY  FORGE  RESEARCH  CENTER. 
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FIGURE  3.4A  IMAGES  OF  HIGH  STREET t  PHOENIXVILLE,  PA,  WITH  AND  WITHOUT 
ADAPTIVE  BEAMFORMING. 
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FIGURE  3.4B  DIFFRACTION-LIMITED  3  cm  RADAR  MAP  OF  SUBURBAN  HOUSING  DEVELOP¬ 
MENT,  PHOENIXVILLE,  PA.,  SHOWING  STREETS  AND  HOUSES.  EACH  PIXEL 
CLUSTER  IS  A  HOUSE.  THE  HOUSES  ARE  ABOUT  20  m  APART.  RESOLUTION 
IS  COMPARABLE  TO  HUMAN  VISION.  THE  ARRAY  LENGTH  IS  83  m. 
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The  second  facility  was  a  physical  aperture  consisting  of  32 
X-band  (3  cm)  receivers  each  with  12  in  horns  (Figures  J  .  5A  and 
B)  .  The  receivers  were  spaced  by  approximately  1  m,  the  element 
placements  being  inaccurate  over  the  32  m  course  of  the  array.  The 
receivers  were  placed  on  a  laboratory  roof  and  were  not  bolted  to 
a  common  platform.  The  upper  portion  of  Figure  3.6  is  a  picture  of 
a  Boeing  727  flying  at  a  distance  of  3.2  km  from  the  laboratory 
made  with  this  array.  The  lower  portion  shows  the  same  data  with¬ 
out  adaptive  beamforming.  The  importance  of  the  self-calibration 
procedure  again  is  evident . 

A  third  facility  is  a  single,  monostatic  X-band  radar  used  in 
an  ISAR  mode.  Figure  3.7  shows  the  pedestal  and  1.2  m  dish 
antenna.  The  flexible  waveguide  connecting  the  antenna  to  the 
equipment  in  the  laboratory  is  visible.  Figure  3 . 7A  is  a  pair  of 
images  made  from  a  common  set  of  data  with  and  without  adaptive 
beamforming.  The  location  of  the  antenna  did  not  change  from  pulse 
to  pulse  in  this  experiment.  The  distortion  instead  was  due  to  the 
perturbations  in  the  flight  path  of  the  airplane  target  due  to 
wind  turbulence  (Figure  3.8).  Adaptive  beamforming  corrected  the 
distortion  (upper  figure) ;  without  it,  the  echo  data  produced  a 
hopeless  jumble  of  pixels. 

3.4  BEAMFORMER  STATISTICS 

Dependence  upon  implanted  corner  reflectors  as  adaptive  beam¬ 
forming  sources  is  unrealistic.  Practical  system  design  requires 
that  echoes  from  targets  themselves  serve  for  the  self-calibration 
process.  Three  studies  have  been  conducted  of  the  frequency  of 
occurrence  of  targets  of  opportunity  that  can  serve  as  beamforming 
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Fig. XT*  32-element  microwave  array.  Each  module  is  a 
microwave  receiver  at  3<m  wavelength.  The  mean  spacing 
is  one  meter.  The  elevation  look  angle  of  each  module  is 
approximately  30°.  Placement  and  pointing  were  done  by 
eye.  This  array  is  highly  distorted  relative  to  the  nominal  0.1 
X  standard  tolerance  in  array  design.  This  distortion  is 
responsible  for  the  poor  imaging  performance  shown  in  Fig.?.  W 
(c);  it  is  electrically  corrected  by  adaptive  beamforming,  as 
is  demonstrated  in  Fig.  (From  [22]  .  ) 


FRONT  VIEW 


REAR  VIEW 


FIGURE  3.5B.  FRONT  AND  REAR  VIEWS  OF  ARRAY 
MODULES.  (FROM  [22].) 
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FIGURE  3.6  MICROWAVE  IMAGE  OF  BOEING  727  OBTAINED 
WITH  32m  X-BAND  RADIO  CAMERA.  AIRPLANE 
'  IS  FLYING  AT  DISTANCE  OF  3.2  Km  FROM 
VALLEY  FORGE  RESEARCH  CENTER.  (From  [22]. 


(a) 

WITH  ADAPTIVE 
BEAMFORMING 


(b) 

PLAN  VIEW  DRAWING 
OF  AIRCRAFT 


(c) 

WITHOUT  ADAPTIVE 
BEAMFORMING 


52 


PHOTO  NOT  AVAILABLE 


FIGURE  3. 


7  MONOSTATIC  X-BAND  RADAR  USED  FOR  ISAR 
IMAGING,  SHOWING  PEDESTAL  AND  1.2m  DISH. 
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FIGURE  3.7A.  MICROWAVE  IMAGE  OF  BOEING  727  OBTAINED  AT  X-BAND 
WITH  MONOSTATIC  RADIO  CAMERA.  AIRPLANE  IS  FLYING 
AT  DISTANCE  OF  2.7  Ko  FROM  VALLEY  FORGE  RESEARCH  CENTER 
(From  [22].) 
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RADAR 


FIGURE  3.8  FLIGHT  PATH  DISTORTION  IN  ISAR. 


sources.  The  first,  by  Yadin,  examined  airborne  radar  data  ob¬ 
tained  by  the  Naval  Research  Laboratory  in  the  upper  UHF  band  (850 
MHz)  [10] .  The  data  were  radar  ground  clutter  obtained  from  an  8- 
element  receiving  system  flown  at  200  knots  over  the  southeastern 
portion  of  the  United  States  at  an  altitude  of  15,000  to  20,000 
ft.  Because  the  antenna  array  was  a  rigid  structure  with  known 
characteristics,  it  was  possible  to  compare  the  adaptively  formed 
radiation  pattern  based  upon  the  ground  clutter  with  the  intrinsic 
pattern  of  the  antenna  itself.  This  comparison  was  made  from 
echoes  from  a  large  number  of  range  bins .  The  criterion  of  whether 
or  not  a  particular  range  bin  contained  a  suitable  adaptive  beam¬ 
forming  source  was  determined  by  the  rise  in  the  sidelobe  level; 
an  acceptable  range  bin  was  one  for  which  the  sidelobe  level  rose 
no  more  than  3  dB.  Based  upon  this  arbitrary  criterion,  it  was 
found  that  approximately  10%  of  the  range  bins  were  acceptable 
candidates  for  adaptive  beamforming.  In  view  of  the  fact  that  at 
most  one  range  bin  per  radar  trace  is  needed  for  the  ABF  process 
and  that  at  the  most  one  trace  per  beamwidth  is  necessary,  the 
observed  frequency  of  occurrence  is  encouraging. 

A  second  study  also  examined  the  frequency  of  occurrence  of 
beamforming  sources  amidst  echoes  from  ground  targets.  The  experi¬ 
ments  are  in  ground-to-ground  imaging  at  very  low  grazing  angles 
[11].  The  array  sketched  in  Figure  3.4  was  employed.  Radiation  was 
in  the  X-band.  Data  were  obtained  at  distances  of  4-7  km  from 
farmland,  and  residential  and  industrial  areas  from  a  town.  The 
results  are  shown  in  Table  3.2.  These  results  are  not  as  opti- 


56 


mistic  as  the  observed  frequency  of  occurrence  in  the  airborne 
case  but  are  nevertheless  very  satisfactory. 


RUN  ,.1XPS 
Industrial 
Residential 
Undeveloped 


BEAMFORMERS 

AVAILABLE 

16.9% 

11 . 6% 

2.7% 


DETECTED 

5.8% 

1.3% 

1.8% 


TABLE  3. 2. VARIANCE  TEST  EFFECTIVENESS  (FROM  [11]).  All  VALUES 

2 

RELATIVE  TO  TOTAL  NUMBER  OF  BINS  PROCESSED.  aNi  <  0.12. 


RANGE  CELL 


1-ZD 

3  m 

%  SUCCESS 

87% 

28% 

NUMBER 

OF 

EXPERIMENTS 

18 

11 

NUMBER 

OF 

SNAPSHOTS 

68 

36 

TABLE  3.3  PERCENTAGE  OF  AIRCRAFT  IMAGING  EXPERIMENTS  IN  WHICH 
ADAPTIVE  BEAMFORMING  WAS  SUCCESSFUL. 
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The  primary  differences  between  the  experiments  were  the 
grazing  angles  and  the  criteria  of  acceptability  of  ABF  range  bin 
data.  The  airborne  grazing  angle  was  a  few  degrees  while  the 
ground-to-ground  grazing  angle  was  nearly  zero.  The  criteria  of 
acceptability  were  also  different.  In  the  airborne  experiment,  a 
range  bin  was  accepted  when  the  radiation  pattern  adaptively 
formed  from  its  echoes  had  satisfactory  sidelobe  properties.  A 
more  sophisticated  measure  was  chosen  for  the  ground-to-ground  ex¬ 
periments.  First,  a  corner  reflector  was  used  to  form  the  weight 
vector  at  the  array,  after  which  an  image  of  a  target  sector  was 
obtained.  This  image  is  the  reference  image.  Second,  the  echoes 
from  each  range  bin  were  used  to  obtain  a  set  of  trial  weight  vec¬ 
tors,  and  an  image  was  formed  from  each  one.  Next,  the  correlation 
coefficient  between  each  trial  image  and  the  reference  image  was 
calculated.  The  criterion  for  acceptability  was  a  correlation  co¬ 
efficient  in  excess  of  0.88,  a  value  obtained  by  a  subjective 
viewing  test  using  63  subjects  [12]. 

Both  tests  lead  to  optimistic  inferences.  Typically  there  are 
hundreds  of  range  bins  in  a  radar  trace.  Thus,  when  the  frequency 
of  occurrence  of  acceptable  range  bins  is  few  per  cent,  the  proba¬ 
bility  is  nearly  unity  that  a  suitable  target  of  opportunity  will 
be  found . 

A  third  study  of  a  different  nature  involves  adaptive  beam- 
forming  and  imaging  on  aircraft,  which  are  isolated  targets.  Here 
the  question  is  the  fraction  of  aircraft  snapshots  obtained  during 
the  passage  of  a  target  through  the  transmitter  beam  for  which  a 
single  range  bin  proves  suitable,  i.e.,  for  which  imaging  can  be 
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performed.  Table  3.3  shows  the  results  of  observations  made  from 
approximately  100  data  sets  obtained  from  about  30  experiments  at 
the  Valley  Forge  Research  Center.  The  results  are  very  satisfac¬ 
tory  for  the  shorter  pulse  duration.  The  reason  that  they  were 
less  satisfactory  for  the  longer  pulse  duration  is  that  the  proba¬ 
bility  that  a  single  dominant  scatterer  exists  in  any  single  range 
bin  decreases  rapidly  with  the  size  of  the  range  cell. 

While  these  results  are  satisfactory  on  the  average,  they  do 
not  guarantee  the  successful  use  of  this  algorithm  at  all  times.  A 
more  robust  procedure  is  described  in  the  next  section  that  in¬ 
creases  the  probability  of  success  still  further.  Moreover,  the 
dominant  scatterer  algorithm  with  the  minimum  amplitude  variance 
test  is  sure  to  fail  under  those  circumstances  in  which  no  domi¬ 
nant  scatterer  exists.  An  example  is  the  backscatter  from  the  sur¬ 
face  of  the  sea.  Fields  of  vegetation,  crowns  of  trees  and  large 
extents  of  sand  are  in  a  similar  category.  In  these  cases,  the 
backscatter  is  statistically  homogeneous  and  a  different  type  of 
algorithm  is  required.  Section  4  addresses  this  problem. 

3.5  MULTIPLE  SCATTERER  ALGORITHM  (MSA) 

When  no  high  quality  adaptive  beamforming  source  is  found,  a 
form  of  range-bin  diversity  sometimes  can  be  employed  to  synthe¬ 
size  a  useful  ABF  source.  The  echo  data  from  several  range  bins 
are  combined  to  produce  a  superior  weight  vector  to  what  would  be 
available  from  any  of  the  individual  range  bins.  A  modification  to 
the  DSA  was  invented  by  Attia  to  deal  with  this  situation  [13]  . 
Experiments  were  reported  in  [14] .  The  method  is  called  by  its  in- 
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ventor  the  robust  minimum  variance  algorithm  (RMVA) ;  here  it  is 
called  the  multiple  scatterer  algorithm  (MSA) . 

Conceptually,  the  procedure  works  as  follows:  in  Figure  3.9 
there  is  a  transmitter  of  beamwidth  A0  illuminating  a  set  of  iso¬ 
lated  scatterers  in  different  range  bins.  These  scatterers  are 
called  candidate  beamformers,  meaning  that  none  is  suitable  for 
adaptive  beamforming  by  itself,  but  that  together  a  satisfactory 
virtual  beamformer  can  be  synthesized.  Amplitude  and  phase  data 
are  recorded  at  each  element  in  the  array  for  each  of  these  range 
bins.  The  algorithm  proceeds  as  follows.  The  range  bin  for  which 

the  echo  amplitudes  at  the  various  receivers  in  the  array  have  the 
smallest  variance  is  selected  as  the  reference  bin  RQ,  as  in  the 

earlier  algorithm.  Second,  the  phases  <|>0^n  of  the  received  sig¬ 
nals  from  RQ  are  subtracted  from  the  phases  of  the  signals  of 

M  candidate  bins,  thereby  removing  the  random  phase  errors  Pn, 
again  as  in  the  former  algorithm.  With  the  |3n  deleted,  the  differ¬ 
ence  phases  form  quadratic  sequences,  modulo  2n,  across  the  array. 

These  phase  sequences  are  then  "unwrapped"  by  removing  the  modulo 
2k  discontinuities.  While  the  unwrapped  phases  are  corrected  for 

the  aperture  distortion,  there  remain  residual  phase  errors  caused 

by  the  clutter.  These  clutter-induced  residual  phase  errors  are 

unbiased  and  therefore  can  be  reduced  by  averaging.  Averaging 
yields  “  n-  0o,n^u^'  which  is  the  averaged  unwrapped  phase  of 

the  difference  phases  of  the  M  candidate  beamformers  from  the  ref¬ 
erence  ABF  bin  at  the  nth  element.  The  superscript  'u'  means  un¬ 
wrapping. 
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ANTENNA  BROADSIDE 


*1  x2  x3  ^ITENNA  ELEMENT  POSITION  *N-1  *N 

FIGURE  3.9  ANTENNA  ARRAY  AND  TARGET  GEOMETRY.  (FROM  [24].) 
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The  steps  can  be  summarized  as  follows: 

1.  Find  R. 

o . 

2.  Form  <j>m,n  -  0o,n  to  remove  (Jn. 

3.  Unwrap  <j>m,n  -  ^0#n  "  <1*0,  n)U  to  remove  modulo  2n 

ambiguities . 

1  M 

4.  Take  the  average  ~  { Z m  ( <J>m  <l>0/n)u  of  the  unwrapped 

phases  to  reduce  the  variance  of  the  clutter-induced 
residual  phase  errors. 

Each  of  the  range  bins  is  assumed  to  contain  a  single  scat- 
terer  embedded  in  clutter  where  the  scatterer  is  not  sufficiently 
large  to  operate  the  DSA.  We  call  such  a  scatterer  subdominant.  In 

none  is  the  target -to-clutter  ratio  adequate  for  the  DSA  to  oper¬ 
ate  satisfactorily.  There  are  M+l  such  scatterers.  Let  Rm,  1  <  m 

<  M  +  l,  be  their  distances  from  the  origin  and  um  =  sin0m  their 
directions  relative  to  the  normal  to  the  array.  The  geometry  is 
shown  in  Figure  3.9.  The  mathematics  of  Section  3.2  can  be  car¬ 
ried  over  directly  from  the  DSA  to  the  MSA  by  replacing  all  sub¬ 
scripts  "o" ,  denoting  range  bin  R0,  by  "m" ,  denoting  range  bin  Rm- 
<t>mf n  is  found  in  this  manner  from  (3.31).  The  difference  phases  are 
obtained  by 

^m,n  “  ^o,n  *  fkxn(um”uo>  “  kxn  *  2Rm  ~  2R0  * 

2  2 
2  um  uo 

+  kxn  *2Rm  ”  2RC  *  +  “  ^o,iJmod  2 it  (3.33) 

which  is  then  unwrapped 
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) 


«fcn,n  ~  <}>o,n)u  =  kxn  (um-uQ) 


+  kx. 


2Rn  *  +  n  “  ^o,n- 


(3.34) 


For  many  microwave  imaging  systems  and  target  geometries  the 
2  2 
2  o  m 

quantity  kx  (:rr—  -  )  «  ft,  even  in  very  near-field  imaging 

n  ZK0  ZKm 

situations.  Further,  1/2RC  =  l/2Rm;  therefore,  (3.34)  is  nearly 

linear  in  x  (first  term)  and  is  perturbed  by  the  distortions  in 
the  phasefronts  caused  by  clutter  in  range  bins  RQ  and  R^  (last  two 

terms)  .  The  average  of  M  such  quantities  will  also  be  a  straight 
line  with  a  smaller  perturbation  provided  that  the  clutter  in  the 
various  candidate  beamformer  range  bins  are  uncorrelated.  The  un¬ 
wrapping  is  illustrated  in  Figure  3.10a  where  the  dots  represent 
^m, n_<^o,n  versus  element  position  n.  Each  of  the  M  curves  <J>m, n  -  <l>o,n 

versus  n  exhibits  discontinuities  at  different  values  of  n. 


Unwrapping  removes  these  discontinuities  and  permits  averaging  the 
M  curves  to  obtain  a  linear  regression  line  with  smaller  perturba¬ 
tions.  The  symbols  at  the  left  and  right  edges  of  Figure  3.10a 

represent  the  unwrapped  values  of  the  symbols  '  .  ’  .  The  residual 
clutter-induced  phase  errors  -  &f>0,n  in  (3.34)  are  depicted  in 

Figure  3.10b. 

The  success  of  the  procedure  depends  upon  the  quality  of  the 

phase  unwrapping.  It  is  clear  that  phase  unwrapping  will  be  a  rel¬ 
atively  easy  task  if  the  slope  of  -  <f>0,n  is  small,  i.e.,  k  (uQ  - 

um)  Ax  is  very  much  less  than  2n  where  Ax  is  the  nominal  interele 
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-  REGRESSION  LINE 


FIGURE  3.10  (a)  PHASE  UNWRAPPING 

(b)  RESIDUAL  PHASE  ERRORS  BETWEEN  THE  REGRESSION  LINE 
AND  THE  UNWRAPPED  PHASES 
(FROM  [24].) 
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ment  spacing.  The  residual  phase  errors  are  easily  obtained  by 
calculating  the  differential  phases  between  the  unwrapped  phases 
and  the  linear  regression  line.  The  dots  in  Figure  3.10b  represent 
5Vn  "  S*o,  n  versus  n.  The  slope  of  a  linear  periodic  array  with 
\/2  spacing  is  7t(um  -  uQ)  .  I  um  -  uQ  |  cannot  exceed  the  transmitter 
beamwidth,  which  is  rarely  more  than  a  few  degrees.  Hence  in  most 

imaging  scenarios  there  will  be  many  data  points  between  the  dis¬ 
continuities  of  <J>m  n  -  <j)Q  n.  However,  for  an  array  with  large  inter¬ 
element  spacing  the  number  of  data  points  between  271  jumps  can  be 
small  and  the  phase  unwrapping  problem  becomes  more  difficult,  in 
which  case  MSA  may  prove  less  applicable.  The  second  requirement 

for  successful  unwrapping  is  that  the  perturbations  in  the 
straight  line  kxn(um  -  uQ)  due  to  clutter  may  be  much  smaller  than 

7t .  It  is  for  this  reason  that  the  range  bin  with  the  smallest 
normalized  amplitude  variance  (and  therefore  the  smallest  clutter- 
induced  phase  perturbation)  is  chosen  to  unwrap  the  phases  of  the 
candidate  beamf ormers . 

The  unwrapped  value  of  <)>mfn  -  <J>Qfn,  denoted  «j>mfn  -  <t>0fn)u,  is 

given  in  (3.34).  The  average  of  M  unwrapped  phase  sequence  curves 
is 

1  M  ^ 

|ln  =  m  X  (<1*111, u  =  kxn  [  (~  X  um)-uQ] 
m=l 


kx 


n 


,  !  X-Um  !-uc 
1  (M  Z  Rm  ]  ~  R0 


+ 


M 


^  5<})mf  n 


o,  n 


(3.35) 
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The  second,  fourth  and  sixth  terms  sum  to  Pn — <}>0  ^  n  (see  (3.31)). 

Adding  (3.31)  to  (3.35)  replaces  these  terms  by  (Jn,  which  reduces 

★ 

the  bias  and  random  errors.  This  sum  is  called  <(>  : 

n 


*n  -  ^n+«>o,n 


kxr 

M 


M 

I 

m=l 


kxn  M  1-um 


+  Pn  +  ~  I  ^n.  (3.36) 


um  - 


M  %  2Rm  rn  '  M  , 

m=l  m  m=l 


The  negative  of  (3.36)  is  the  MSA  error-compensating  phase  shift 
introduced  into  the  nth  channel.  Therefore,  the  nth  component  of 
the  MSA  weight  vector  is 

wn,  msa  =  exp(-j<(»n) 


(3.37) 


The  first  two  terms  in  (3.36)  point  and  focus  the  array  to 
map  coordinates  (u,R) : 


l  M  i 

{M  ^  Um/  M 

m=l 


M 


I 

m=l 


The  third  term  is  the  phase  error  caused  by  electrical  phase  dis¬ 
tortion  and  element  position  error.  The  fourth  term  is  the  aver¬ 
aged  clutter-induced  pha^e  error  of  the  M  candidate  beamformers. 
This  term  is  smaller  than  the  corresponding  DSA  term  (8<j>0  n  in 

(3.31))  because  of  the  averaging. 

3.6  COMPARISON  OF  DSA  AND  MSA 

In  this  section,  the  angular  biases  and  errors  in  the  DSA 

and  MSA  are  compared.  Equations  (3.31  and  (3.36)  define  the  error- 

compensating  phase  shifters.  The  first  term  in  (3.31)  is  kxnuQ  and 

M 

in  (3.36)  it  is  (kxn£um) /M.  These  terms  steer  the  beams.  The  DSA 

m 

centers  the  phase-corrected  image  on  uQ  whereas  the  MSA  centers  the 
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phase-corrected  image  on  (Emum) /M,  which  is  the  average  of  the 

locations  of  the  M  candidate  beamformers.  The  average  is  more 
likely  to  be  closer  to  the  center  of  the  FOV  than  the  location  of 
the  reference  beamformer  uD.  Therefore,  one  advantage  of  the  MSA 
over  the  DSA  is  that  it  focuses  the  array  closer  to  the  center  of 
the  FOV  than  the  DSA,  thereby  reducing  the  uncertainty  in  target 
location . 

An  estimate  of  the  improvement  in  angular  accuracy  is  ob¬ 
tained  using  the  theorem  in  statistics  that  the  standard  deviation 
of  the  sample  mean  of  M  independent  random  variables  varies  as 
M1/2,  and  the  realistic  assumption  that  the  angular  locations  of 

the  candidate  beamformers  are  random  and  independent .  The  expected 
location  of  uQ  is  within  the  main  lobe  (of  width  A0)  of  the  trans¬ 
mitting  beam  for  otherwise  the  target  would  be  too  weakly  illumi¬ 
nated  to  become  the  reference  reflector.  For  the  same  reason  the 
directions  to  all  the  other  candidate  beamformers  also  fall  within 
the  main  beam.  These  bearings  to  the  targets  may  be  assumed  to  be 
approximately  uniformally  distributed  within  the  main  beam  because 
there  is  no  physical  reason  for  their  distribution  to  be  biassed. 
Thus,  after  averaging,  the  standard  deviation  of  the  sample-mean 
pointing  direction  is  Gq/M1^2.  Hence  the  accuracy  in  the  angular 

coordinate  of  the  microwave  image  is  improved  by  the  factor  M1/2. 

The  second  terms  in  (3.31)  and  (3.36)  focus  the  array  to  Rc 
and  (~  Z  R^)-1/  respectively.  There  is  no  advantage  of  one  over 

the  other  because  the  quadratic  phase  curvature  from  either  dis¬ 
tance  is  deterministic  and  therefore  can  be  removed  to  permit  re- 
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focusing  the  array  for  range  bins  lying  outside  the  depth  of  field 


of  the  array. 

The  third  terms  in  (3.31)  and  (3.36)  are  the  phase  errors 
caused  by  electrical  phase  distortion  and  element  position  error. 
This  error  is  constant  for  all  range  bins  and  is  corrected  by  both 
ABF  procedures . 

The  fourth  terms  in  (3.31)  and  (3.36)  are  the  clutter-induced 

phase  errors.  The  MSA  term  is  smaller  because  of  the  M-fold  aver¬ 
aging.  Let  Odsa  be  the  standard  deviation  of  8<|>0  n  of  (3.31)  and 

l  M 

aMSA  be  the  standard  deviation  of  ~  £  &|>m, n  of  (3.36)  .  Given  the 

m=l 

realistic  assumption  that  all  phase  errors  are  independent,  and 
under  the  assumption  that  the  S^n  and  8<J>0fn  have  common  variances, 

the  relation  between  the  phase-error  variances  is 


The 

radiation 


°MSA  *  M  °dsA 

significance  of  (3.38) 

2  ^  2 

pattern.  Omsa  and  <*DSA 


(3.38) 

lies  in  its  influence  upon  the 
are  the  variances  of  the  phase 


errors  in  the  phase  corrections  across  the  array  after  ABF.  It  is 
these  residual  errors  that  determine  the  qualities  of  the 
radiation  patterns  formed  by  these  algorithms.  It  has  been  shown 

that  the  expected  power  pattern  Eff*  when  phase  errors  are  present 

and  amplitude  errors  are  not,  normalized  so  that  the  on-axis  (u=o) 

gain  in  the  absence  of  phase  errors  is  unity,  is  [after  (13.7)  of 


[9]  ) 


Eff"  *  e-°2  f0f0*  +  —  TT 


I?  o2 


a2 


N 


For  the  problem  being  discussed  C2  is  either  0DgA  or  a, 


(3.39) 
2 

■’msa  ' 
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f  =  f (u)  is  the  radiation  pattern  of  the  array  with  phase  error, 

I? 

fc  =  fc(u)  is  the  design  pattern,  —  is  the  square  of  the  ratio  of 

32 

the  rms  to  mean  value  of  the  aperture  weights  and  is  normally 
close  to  and  slightly  greater  than  unity,  and  N  is  the  number  of 
elements.  Clearly  O2  >  0  reduces  main  lobe  gain  and  increases 
sidelobes.  Only  when  <JDSA  is  not  sufficiently  small  is  there  reason 

to  reduce  the  residual  phase-error  variance  with  the  more 

2  , 

complicated  algorithm.  For  example,  assume  that  gdsa  =  0.1  (rad) 


N  =  300,  the  allowed  loss  in  main  lobe  gain  is  0.5  dB  and  the 
sidelobe  specification,  after  ABF,  is  -30  dB.  The  expected  main 
lobe  gain,  relative  to  the  design  gain,  is  e*0-1  *  0.9  or  -0.43  dB, 

which  is  acceptable.  However,  the  sidelobe  level  is  too  high.  The 
average  sidelobe  is  approximately  0.1/300  or  about  -35  dB  and  the 
peak  sidelobe  is  typically  10  dB  higher,  which  exceeds  the 
specification  by  5  dB.  Equation  (3.38)  indicates  that  MSA  will 
meet  the  specifications  for  a  value  of  M  =  3. 

However,  (3.38)  is  somewhat  optimistic  because  the  second 
assumption,  equality  of  variances,  is  not  often  satisfied.  The 
"best"  range  bin  is  selected  for  the  DSA,  implying  that  the  rms  of 
the  phase  perturbations  54>mfn  of  the  MSA  exceeds  that  of  the  5<)>o,n 
of  the  DSA.  Let  a  >  1  be  the  average  of  the  variances  in  the  M 
candidate  beamforming  bins  relative  to  the  DSA  reference  ABF  bin. 
Then  the  proper  relation  between  the  phase-error  variances  of  the 
two  algorithms  becomes 
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(3.40) 


2  _  a  2 

aMSA  “  M  °DSA 


Because  MSA  is  a  more  complicated  algorithm,  achievement  of 
the  improvement  given  by  (3.40)  requires  additional  computation. 
Experience  with  the  algorithm  indicates  that  the  computation  is  2- 
3  times  larger  than  the  DSA  process. 

MSA  plays  the  role  of  a  threshold  extender  to  DSA.  The  lowest 
amplitude  variance  bin  is  selected  for  the  DSA  reference.  Only  if 
it  proves  inadequate  is  the  more  complicated  MSA  employed.  The 
reference  bin  data  are  used  to  unwrap  the  phases  of  a  small  set  of 
the  M  candidate  beamformers.  A  few  dB  reduction  in  the  phase-error 
variance  has  been  observed  in  ground  clutter  experiments.  This 
advantage  is  demonstrated  in  the  next  section. 

The  DSA  and  MSA  can  also  be  compared  by  correlating  the  DSA 
and  MSA  phase-corrected  images  with  a  perfect  image.  The  correla¬ 
tion  measure  has  proved  useful  in  comparing  microwave  images  [12]. 

The  correlation  between  two  images  I^u)  and  I2 (u)  is 

1 

j  Ii(u)  I2 (u) du 

P"  I  2  i  2  (3'41) 

[  J  I1(u)du  J  I2(u)du]1/2 

-1  -1 

where  (u)  is  the  error-free  image  of  the  target  scene  and  l2(u) 

is  the  erroneous  image.  Assume  that  the  perturbations  in  the  phase 
corrections  caused  by  clutter  are  independent  and  identically  dis¬ 
tributed  Gaussian  random  variables  with  zero  mean  and  variance  a2. 


A  lower  bound  on  the  correlation  p between  the  error-free  and  erro- 
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neous  images  has  been  derived  in  [16]  in  terms  of  the  residual 

phase-error  variance;  it  is  given  by 

po=  j  (1  +  e_°2  )  =  e"°2/2  (3.42) 

so  that  p  >  po  By  comparing  (3.42)  and  (3.39)  we  see  that  the 
expected  value  of  the  mainbeam  gain  relative  to  the  design  gain  GQ 
is  the  square  of  po: 

E  =  exp("V  *  Po  <3-43> 

3.7  DSA  AND  MSA  EXPERIMENTS 

The  83-m,  330-element  quasi-linear  array  pictured  in  Figure 
3.4  was  used  to  image  the  nuclear  power  plant  in  Limerick,  PA, 
located  17.6  km  from  the  array.  The  array  was  highly  distorted; 
the  locations  of  the  receivers  differed  from  their  assumed  loca¬ 
tions  by  one  or  more  wavelengths  in  each  dimension. 

Beamforming  bins  having  normalized  amplitude  variances  not 
larger  than  0.12  have  been  shown  to  yield  satisfactory  images  when 
the  DSA  is  used  for  phase  cohering  an  array  [11] .  A  survey  of  the 
normalized  amplitude  variances  of  the  75  range  bins  in  the  FOV 
yielded  6  bins  that  could  possibly  contain  a  single  subdominant 
reflector  embedded  in  clutter.  These  range  bins  were  ordered  in 
accordance  with  their  variances  and  numbered  0  through  5;  their 
normalized  amplitude  variances  were  0.06,  0.07,  0.09,  0.10,  0.10 

and  0.12.  Bin  0  contains  the  highest  quality  beamformer.  A  compar¬ 
ison  was  made  between  the  DSA  and  the  MSA  both  using  bin  0  as  the 
reference  bin.  Bins  1  and  2,  the  next  highest  qualitiy  bins,  were 
used  as  the  candidate  bins  for  the  MSA  test.  These  results,  shown 
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compare  the  DSA  and  the 


in  parts  a  and  b  of  Figures  3.11-3.13, 

MSA  under  low  variance  conditions,  i.e.,  highly  favorable  for 
imaging . 

A  second  comparison  was  made  using  high  variance  range  bins, 
none  of  which  was  favorable  for  imaging.  These  range  bins  had  nor¬ 
malized  amplitude  variances  of  0.20,  0.21  and  0.23  and  were  num¬ 
bered  6  through  8.  The  results  of  these  experiments  are  shown  in 
parts  c  and  d  of  Figures  3.11-3.13. 

3.7.1  1-D  CROSS-RANGE  IMAGES  OF  BIN  2 

Figure  3.11  shows  the  cross-range  images  of  bin  2  obtained 

from  the  four  experiments.  The  dotted  line,  which  is  25  dB  below 

the  peak,  is  a  typical  threshold  value  below  which  signals  are  not 

displayed  in  range-azimuth  imagery.  It  is  close  to  the  average 
sidelobe  level  of  10  log  where  N  is  the  number  of  antenna  array 

elements.  The  upper  two  parts  compare  the  low  variance  experiments 
and  the  lower  two  parts  compare  the  high  variance  experiments. 

Figure  3.11a  shows  that  the  DSA  peak  sidelobe  level  (PSL)  is 
about  -13  dB.  Figure  3.11b  shows  the  improvement  due  to  the  MSA: 
the  PSL  is  decreased  to  about  -16  dB.  Furthermore,  the  sidelobe 
artifacts  around  the  main  beam,  evident  in  the  DSA  image,  have 
been  removed. 

The  high  variance  comparison  is  made  in  the  lower  two  fig¬ 
ures.  Figure  3.11c  shows  the  performance  of  the  DSA  with  bin  8  as 
the  reference  bin.  The  PSL  has  risen  to  about  -7  dB  and  many  side- 
lobes  rise  above  the  threshold.  This  experiment  demonstrates  that 

bin  8,  whose  normalized  amplitude  variance  is  0.23,  is  a  very  poor 
beamformer  for  the  DSA.  However,  when  additional  candidate  beam- 
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FIGURE  3.11  CROSS-RANGE  IMAGE  OF  BIN  2  OF  LIMERICK  NUCLEAR  POWER  PLANT 
ICW  VARIANCE  CASE 

(a)  DSA  WITH  BIN  0  AS  THE  REFERENCE  BIN 

(t)  MSA  WITH  BIN  0  AS  THE  REFERENCE  BIN  AND  BINS  1-2  AS  THE  CANDIDATE  BINS 

HIGH  VARIANCE  CASE 

(c)  DSA  WITH  BIN  8  AS  THE  REFERENCE  BIN 

(2)  MSA  WITH  BIN  8  AS  .TOE  REFERENCE  BIN  AND  BINS  -6-7  AS  THE  CANDIDATE  BINS 

(FROM  [24].) 


former  bins  of  equally  poor  quality  are  added,  the  MSA  performs 
reasonably  satisfactorily.  Figure  3. lid  shows  this  case.  The  PSL 
has  decreased  to  about  -11  dB.  This  MSA  image  is  almost  as  good  as 
the  DSA  image  of  Figure  3.11a  formed  from  the  highest  quality 
beamformer . * 

In  addition  to  a  comparison  based  upon  the  PSL,  a  second  use¬ 
ful  comparison  is  the  number  of  apparent  targets  appearing  above 
the  display  threshold.  Both  values  are  listed  in  Table  3.4  for  the 
four  cases,  along  with  the  normalized  amplitude  variance  informa¬ 
tion.  From  the  table  it  is  evident  that  the  MSA  high  variance  case 
approaches  the  quality  of  the  low  variance  DSA  case. 

Figure  3.11  shows  the  details  of  the  residual  phase  errors 
across  the  aperture  for  the  four  cases.  The  rms  errors  and  the 
losses  in  mainbeam  gain  (calculated  from  (3.43))  are  listed  in 
Table  3.5. 

3.7.2  2— D  RANGE-AZIMUTH  IMAGES 

Figure  3.13  shows  range-azimuth  images  of  buildings  at  the 

site  of  the  Limerick  nuclear  power  plant.  The  upper  two  images, 
both  obtained  from  low  variance  beamformers,  are  very  similar  and 
both  are  satisfactory.  Although  the  MSA  reduces  the  phase  errors 
and  thereby  intrinsically  produces  a  higher  quality  image,  the  DSA 
image  of  Figure  3.13a  is  of  sufficient  quality  so  that  use  of  the 
MSA  is  not  warranted. 

♦It  is  interesting  to  observe  the  shift  in  the  optical  axes  of  the 
lower  images  relative  to  the  upper  ones.  The  radio  camera  points 
in  the  direction  of  the  reference  reflector.  The  shift  reflects 
the  difference  in  directions  of  the  reference  reflectors  in  bins  0 
and  8 . 
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FIGURE  3.12.  PHASE  ERRORS  OF  BIN  2  OF  LIMERICK  NUCLEAR  POWER  PLANT 


LOW  VARIANCE  CASE 

(a)  DSA  WITH  BIN  0  AS  TOE  REFERENCE  BIN 

Cb)  MSA  WITH  BIN  0  AS  TOE  REFERENCE  BLN  AND  BINS  1-2  AS  THE  CANDIDATE  BINS 
HIGH  VARIANCE  CASE 

(c)  DSA  WITH  BIN  8  AS  TOE  REFERENCE  BIN 

(d)  MSA  WITH  BIN  8  AS  TOE  REFERENCE  BLN  AND  BINS  6-7  AS  THE  CANDIDATE  BINS 
(FROM  [24].) 
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NOV.  87 


VFRC-1 7-87 
meters 


FIGURE  3.13  RANGE- AZIMUTH  IMAGES  AT  SITE  OF  LIMERICK  NUCLEAR  POWER  PLANT. 

LOW  VARIANCE  CASE 

(a)  DSA  WITH  BIN  0  AS  THE  REFERENCE  BIN 

(b)  MSA  WITH  BIN  0  AS  THE  REFERENCE  BIN  AND  BINS  1-2  AS  THE  CANDIDATE  BINS 

HIGH  VARIANCE  CASE 

(c)  DSA  WITH  BIN  8  AS  THE  REFERENCE  BIN 

(d)  MSA  WITH  BIN  8  AS  THE  REFERENCE  BIN  AND  BINS  6-7  AS  THE  CANDIDATE  BINS 

(FROM  [24].) 
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LOW  VARIANCE 

HIGH 

VARIANCE 

SELF -CAL I BRAT I ON  ALGORITHM 

DSA 

MSA 

DSA 

MSA 

EXPERIMENT  IN  FIGURE  3.11 

a 

b 

d 

d 

AMPLITUDE  VARIANCE 

REFERENCE  BEAMFORMER 

0.06 

0.06 

0.23 

0.23 

CANDIDATE  BEAMFORMERS 

0.07 

0.20 

0.09 

0.21 

PEAK  SIDE  LOBE  (dB) 

-13 

-16 

-7 

-11 

NUMBER  OF  TARGETS  >  -25  dB 

15 

10 

26 

15 

TABLE  3.4  COMPARISON  OF  DSA  AND 

MSA  SELF-CALIBRATING 

ALGORITHMS 

IN  LOW  AND 

HIGH  VARIANCE  CASES.  THE  MSA  HIGH  VARIANCE  CASE  IS  ALMOST  AS  GOOD 
AS  THE  DSA  LOW  VARIANCE  CASE. 


LOW  VARIANCE  HIGH.  VARIANCE 


SELF-CALIBRATION  ALGORITHM 

DSA 

MSA 

DSA 

MSA 

EXPERIMENT  IN  FIGURE  3.12 

a 

b 

c 

d 

RESIDUAL  RMS  PHASE  ERROR  (deg) 

23.57 

10.07 

46.74 

33.35 

AG  (dB) 

-0.73 

-0.13 

-2.89 

-1.47 

TABLE  3.5.  RESIDUAL  RMS  PHASE 

ERRORS  FOR 

THE  FOUR  CASES 

OF  FIGURE 

3.11  AND 

THE  LOSSES  IN  MAIN  BEAM  GAIN  THAT  RESULT  FROM  THESE  ERRORS. 
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Not  so  with  the  imagery  from  the  low  quality  beamformers. 
Figure  3.13c  shows  the  DSA  image  driven  by  a  range  bin  with  a  nor¬ 
malized  amplitude  variance  of  0.23.  This  image  is  badly  degraded 
by  high  sidelobes;  it  is  not  a  satisfactory  image.  The  same  low 
quality  reference  beamformer  plus  two  equally  low  quality  candi¬ 
date  bins  produce  a  considerably  improved  MSA  image  (Figure 
3.13d).  Here  the  sidelobe  artifacts  have  been  largely  removed  and 
the  image  is  not  significantly  inferior  to  the  high  quality  images 
in  Figures  3.13a  and  3.13b. 

3.7.3  IMAGE  CORRELATION 

It  is  interesting  to  form  an  image  correlation  comparison 
based  on  the  definition  given  in  (3.41).  The  best  of  the  four  im¬ 
ages  is  Figure  3.13b,  as  determined  by  the  residual  RMS  phase 
data.  Using  it  as  the  reference  image,  its  correlations  with 
Figures  3.13a,  c  and  d  are  0.98,  0.87  and  0.92,  respectively.  In 
[12]  it  was  shown  that  a  correlation  coefficient  of  0.9  or  greater 
indicates  satisfactory  image  formation.  Based  on  that  criterion, 
experiment  c  fails  while  a  and  d  are  satisfactory.  Thus  we  again 
see  that  three  range  bins  having  normalized  amplitude  variances 
0.20  or  greater  can  be  used  with  the  MSA  to  obtain  a  useful  image, 
one  nearly  as  good  as  the  image  obtained  using  the  DSA  with  the 
best  beamformer.  Moreover,  if  only  one  of  the  three  range  bins  is 
used  for  beamforming  purposes,  the  resulting  image  will  have  large 
sidelobe  artifacts,  as  is  demonstrated  in  Figures  3.11c  and  3.13c. 
In  addition,  the  loss  in  mainlobe  gain  relative  to  experiment  b 
can  be  calculated  using  (3.43).  Table  3.6  summarizes  these  find¬ 
ings  . 
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LOW  VARIANCE 


HIGH  VARIANCE 


SELF-CALIBRATION  ALGORITHM 

DSA 

MSA 

DSA 

MSA 

EXPERIMENT  IN  FIGURE  3.13 

a 

b (REFERENCE) 

C 

d 

IMAGE  CORRELATION 

0.98 

1.0 

0.87 

0.92 

AG(dB) 

-0.18 

0 

-1.21 

-0.72 

TABLE  3.6  IMAGE  CORRELATIONS  AND  LOSSES  OF  MAIN  BEAM  GAIN  RELATIVE  TO 

EXPERIMENT  b  AGAIN  DEMONSTRATING  THAT  THE  MSA  HIGH  VARIANCE  CASE 
(FIGURE  3.13d)  IS  NEARLY  AS  GOOD  AS  THE  DSA  LOW  VARIANCE  CASE 
(FIGURE  3.13a) . 

The  last  experiment  illustrates  the  capability  of  a  dominant 
scatterer  algorithm  to  image  a  distant  target  with  human-optical 
resolution.  The  VFRC  83-m  array  was  used  to  image  an  industrial 
site  in  Phoenixville,  PA,  located  at  a  distance  of  8.2  km  from  the 
array.  Figure  3.14a  is  an  aerial  photo  of  the  Cromby  Power  Plant 
which  shows  two  high  stacks,  a  water  tower,  the  power  plant  build¬ 
ing,  and  two  connected  storage  silos.  The  microwave  image  of  the 
white  blocked  area  in  the  upper  figure  is  shown  below.  The  radar 
line  of  sight  is  indicated  by  the  white  vertical  line  in  the  cen¬ 
ter  of  the  upper  figure.  The  radio  camera  image  was  obtained  by 
applying  the  MSA  to  the  best  set  of  three  range  bins;  it  reveals 
the  two  high  stacks  (A  and  B) ,  the  water  tower  (C)  ,  the  two  con¬ 
nected  silos  (D) ,  and  the  power  plant  building  (E) .  The  resolution 
cell  is  3  m  in  both  range  and  cross  range.  The  MSA  phase  correc¬ 
tions  compensated  for  the  large  geometric  distortion  of  *.he  phased 
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array  antenna,  and  very  little  sidelobe  energy  is  evident  in  the 
microwave  image . 

The  variance  range  from  0.12  to  0.20  is  the  DSA  transition 
zone.  The  algorithm  always  works  at  the  low  end  and  never  works  at 
the  high  end.  The  MSA,  on  the  other  hand,  provides  a  useful  image 
from  three  range  bins  having  normalized  amplitude  variances  as 
large  as  the  order  of  0.2.  In  this  sense  MSA  acts  like  a  threshold 
extender,  permitting  the  radio  camera  to  function  under  adverse 
clutter  conditions. 

The  performance  of  each  algorithm  can  be  improved  by  subarray 
processing . 

3.8  USE  OF  SUBARRAYS 

Other  variations  of  the  basic  algorithm  improve  its  opera¬ 
tion.  One  modification  is  the  division  of  the  large  array  into 
sections  or  subarrays  to  enhance  the  adaptive  beamforming  process 
when  the  beamforming  source  is  not  an  ideal  point  reflector. 
Figure  3.15  shows  the  problem.  A  target  of  size  T  reradiates  ener¬ 
gy  toward  the  array.  The  radiation  is  always  lobular  and  the  nomi¬ 
nal  lobe  width  is  X/ T  almost  independent  of  the  distribution  of 

the  scatterers  on  the  target.  The  phase  across  the  central  portion 
of  a  lobe  is  fairly  constant,  while  it  changes  abruptly  and  often 
randomly  from  lobe  to  lobe.  For  the  ABF  process  to  proceed 
smoothly  it  is  necessary  for  the  entire  array  to  be  illuminated  by 
the  central  portion  of  a  single  lobe.  A  typical  lobe  cross  section 
at  the  array,  which  is  a  distance  R  from  the  target,  is  Xr/T. 
Unless  L,  the  array  size,  is  smaller  than  Xr/T,  this  condition  is 
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FIGURE  3".15  The  reradiation  pattern  of  the  target 
is  lobular.  (From  [4].) 


not  likely  to  be  satisfied.  It  is  shown  in  [4]  and  [S1  that  the 
minimum  target  size  for  satisfactory  beamforming  is 

Xr 

T  <  tj^jT  (3.44) 

Experience  with  manmade  reflectors  indicates  that  the 
condition  (3.44)  is  violated  sometimes  and  imaging  quality  with 
the  DSA  or  the  MSA  degrades.  It  is  often  possible  to  correct  this 
situation  by  breaking  the  array  into  sections  or  subarrays  each 
small  enough  to  satisfy  (3.44) .  Let  <  L  be  such  a  subarray  such 

that 

Xr 

T  <  ^  (3.45) 

The  phase  of  the  reradiated  wavefront  from  the  target  across  the 
subarray  is  nearly  constant  and  phase  corrections  for  element 
position  errors,  electrical  mistunings,  etc.,  can  be  accomplished. 
Let  L2 ,  L 3,  etc.,  be  other  such  subarrays.  Phase  compensation  is 

performed  across  each  subarray.  Thus  each  subarray  has  its  aberra¬ 
tions  corrected.  However,  the  phase  reference  for  a  particular 
subarray  is  the  phase  of  the  lobe  of  the  reradiation  pattern  of 
the  target  that  illuminates  that  subarray.  Because  the  lobe  phases 
must  be  assumed  to  be  random  with  respect  to  each  other,  the  array 
segments,  each  internally  phase  synchronized,  are  also  randomly 
phased  with  respect  to  each  other.  Thus  the  remaining  problem  is 
to  link  the  array  segments  by  phase  correcting  for  the  random 
phase  differences  between  segments  of  the  array.  Figure  3.16 
illustrates  this  situation.  The  abscissa  is  range  and  the  ordinate 
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RANGE 


FIGURE  3.16  SUBARRAYS  L^,  Lg  AND  Lc  SELF-COHERE  ON 

TARGETS  A,  B  AND  C.  LA  AND  Lg  ARE  OVERLAPPING 
SUBARRAYS;  Lg  AND  Lc  ARE  NONOVERLAPPING. 
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is  element  number.  Target  echoes  fall  in  the  shaded  area.  Let  A,  B 

and  C  be  prominent  scatterers.  None  persists  throughout  the  array, 
however.  A  provides  a  strong,  stable  echo  over  that  portion  of 

the  array  labeled  LA,  and  similarly  for  B  and  C.  Subarrays  LA  and 

Lb  overlap;  LB  and  Lc  do  not. 

Array  elements  within  LA  are  self-calibrated  by  the  echoes 

from  A  using,  say,  the  DSA.  Echoes  from  B  self-calibrate  subarray 
Lb.  The  self-calibrated  subarrays  do  not  share  a  common  phase  ref¬ 
erence  and  therefore  exhibit  a  random  offset  phase  between  them. 
The  offset  between  LA  and  LB  can  easily  be  measured  because  targets 

A  and  B  appear  in  both  data  sets.  Let  be  the  phase  of  target  A 
as  observed  by  subarray  LA  and  <|>AB  its  phase  as  observed  by  sub¬ 
array  Lb.  The  difference  (j)^  -  <|>AB  is  the  phase  offset  between  the 

self-calibrated  subarrays.  Adding  this  difference  to  all  the 
phases  in  the  LB  data  set  eliminates  the  random  phase  offset  and 

extends  the  calibration  over  the  two  subarrays. 

An  additional  step  is  required  when  the  subarrays  do  not 
overlap.  In  that  case  a  target  common  to  both  subarrays  is  found 
after  self-calibration.  Call  this  target  D.  The  phase  difference 
<t>E>B~  0DC  Is  added  to  all  the  phases  in  the  Lc  data  set  to  eliminate 

the  random  phase  offset  between  subarrays  B  and  C. 

The  technique  is  more  complex  when  applied  to  the  MSA.  Two 
useful  procedures  for  solving  the  problem  are  now  discussed.  The 
first  deals  with  non-overlapping  subarrays,  the  second  with  over¬ 
lapping  array  segments.  Both  techniques  are  described  in  the  con¬ 
text  of  the  MSA. 

The  phase  sequence  <j>0>n  from  the  reference  range  bin  RQ^  given 


in  (3.31),  is 
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♦  Sfo.n  +  Pn 


(3.46) 


^o,  n 


kxnuo 


2RC 


where  5<j>0fn  are  the  perturbations  in  the  phasefront  caused  by 
clutter  and  |3n  are  the  phase  errors  due  to  errors  in  array 

geometry  and  mistunings  in  the  receivers.  The  averaged  unwrapped 
phase  differences  of  the  sequences  <j)m^  n  from  M  candidate 

beamformers,  given  in  (3.35),  are 


1  M 

—  —  51  (^,11  -$o,n)u 

M 

m=l 

.  2  .  2  _  2 

kx  1-un 

kxn(-  I  um-uQ)  2  (M  2  Rm  R0  ) 

+  2  8^,11)  “  (3.47) 


The  nth  component  of  the  weight  vector  formed  by  the  MSA  is  (from 
(3.37)  ) 


wn,MSA  -  exp[-j  <m  +  80,n>  1  (3.48) 

Eq.  (3.47)  is  nearly  a  straight  line  (because  the  quadratic 
term  is  very  small)  perturbed  by  zero-mean  clutter-induced  phase 
errors.  Such  a  sequence  is  illustrated  in  Figure  3.10a.  This 
figure  portrays  the  case  in  which  the  main  reflector  is  physically 
small,  i.e.,  satisfies  condition  (3.44),  and  overwhelmingly 
dominates  the  backscatter  from  the  clutter.  In  this  case  the  data 
points  cluster  about  a  linear  regression  line. 

Sometimes,  however,  condition  (3.44)  is  not  satisfied.  Then, 
due  to  the  lobular  nature  of  the  reradiation,  the  data  sequence 
exhibits  the  phenomenon  shown  in  Figure  3.16.  The  low  spatial 
frequency  variation  shown  there  induces  a  phase  error  that  adds  to 
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an  independent,  random  perturbation  from  element  to  element.  This 
phenomenon  is  evident  in  the  residual  phase  errors  shown  in  Figure 
3.17  and  indicates  a  lobular  reradiation  across  the  array.  The 
residual  phase-error  variance  is  increased  because  it  is  the  sum 
of  the  variances  of  the  random  component  and  the  low  frequency 
undulation.  It  is  often  possible  to  remove  much  of  the  latter 
contribution  by  estimating,  from  the  data,  the  phase  change  from 
lobe  to  lobe.  These  phase  jumps  are  then  deleted  from  the  data. 

Figure  3.17  illustrates  one  procedure.  A  jin  sequence  is 

illustrated  in  the  upper  figure.  The  array  is  broken  into  three 
subarrays  and  a  linear  regression  line  (least  mean  square  fit  to 
the  data)  is  found  for  each  subarray  (a)  .  The  result  is  a  lower 
residual  phase  error  variance  (b)  .  The  discontinuities  in  the 
piecewise  linear  fit  of  (a)  to  the  data  are  assumed  to  be  due  to 
phase  changes  between  lobes  in  the  reradiation  field;  they  are 
removed  by  adding  fixed  phases  to  each  section  to  make  the 

piecewise  linear  fit  continuous  (c)  .  The  modified  sequence  of 

* 

averaged  unwrapped  differential  phases  is  called  p.  and  the  nth 

n 

component  of  the  modified  MSA  weight  becomes 

“n,MSA  “  exp  Z  ['j(pn  +  ♦».■>>]  <3-49> 
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FIGURE  3.17  NONOVERLAPPING  SUBARRAY  PROCESSING. 

(a)  REGRESSION  LINE  IN  EACH  SUBARRAY. 

(b)  RESIDUAL  PHASE  ERRORS  AFTER  REMOVING  THE 
REGRESSION  LINE  FROM  THE  DIFFERENTIAL  PHASES. 

(c)  VERTICALLY  CORRECTED  DIFFERENTIAL  PHASES  FROM  (a). 
(FROM  (24]  AND  APPENDIX.) 
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Figure  3.18  illustrates  another  subarray  procedure  for  esti¬ 
mating  the  correction  phases.  The  arrays  are  overlapped  in  this 
case.  Exactly  as  described  before,  regression  lines  are  fit  to 
the  subarray  Hn  data  (a) .  Each  overlap  region  now  has  two  regres¬ 
sion  lines.  Each  pair  is  averaged  to  form  a  new  regression  line  in 
its  limited  interval.  The  result  is  a  discontinuous  piecewise  lin¬ 
ear  fit  to  the  data,  which  are  then  adjusted  exactly  as  described 

before  to  obtain  a  continuous  piecewise  linear  fit  (b) .  The  phase- 

* 

corrected  data  are  the  \l  and  the  modified  MSA  weight  vector  is 

n 

given  by  (3.49). 

3.9  SUBARRAY  EXPERIMENTS 

Figure  3.19  shows  four  1-D  images  of  a  corner  reflector  for 
which  beamforming  was  accomplished  with  a  relatively  unsatisfac¬ 
tory  scatterer,  identified  in  the  figure  as  the  echo  from  bin  61. 
The  target  area  was  a  residential  street  in  Phoenixville, 
Pennsylvania.  The  amplitude  variance  from  that  range  bin  was  0.21. 
Figure  3.19a  shows  the  DSA  performance.  The  maximum  amplitude 
reached  is  53.26  units  and  the  peak  sidelobe  is  -6.9  dB.  The  image 
extends  over  nearly  30  mrad  at  the  -25  dB  level  and  does  not  look 
as  expected  from  a  corner  reflector. 

Figures  3.19b,  c  and  d  show  the  performances  of  MSA,  MSA  with 
non-overlapping  subarrays  and  MSA  with  overlapping  subarrays, 
respectively.  The  superiority  of  Figure  3.19d  is  evident.  The  peak 
amplitude  has  increased  to  69.45  units,  indicating  an  increase  of 
2.3  dB  in  main  lobe  gain,  and  the  peak  sidelobe  has  dropped  to  - 
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FIGURE  3.18  OVERLAPPING  SUBARRAY  PROCESSING. 

(a)  REGRESSION  LINE  IN  EACH  SUBARRAY. 

(b)  RESIDUAL  PHASE  ERRORS  AFTER  REMOVING  THE 
REGRESSION  LINE  FROM  THE  DIFFERENTIAL  PHASES. 

(c)  CORRECTED  RESIDUAL  PHASE  ERRORS. 

(d)  CORRECTED  REGRESSION  LINE  AND  CORRECTED  PHASES. 
(FROM  124]  AND  APPENDIX.) 
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MSA-NSP 


FIGURE  3.19  CROSS-RANGE  IMAGES  OF  CORNER  REFLECTOR.  (FROM  [24]  AND  APPENDIX.) 

(a)  DSA  WITH  BIN  61 

(b)  MSA  WITH  BINS  61,  67  and  68 

(c)  MSA  WITH  NONOVERLAPPING  SUBARRAYS,  USING  BINS  61,  67,  68  AND 
THREE  SUBARRAYS 

(d)  MSA  WITH  OVERLAPPING  SUBARRAYS,  USING  BINS  61,  67,  68  AND 
THREE  SUBARRAYS 


14.4 

dB.  The  extent  of 

the  image 

above  the 

-25  dB  level  is  reduced 

to  17 

mrad.  Table  3.7  lists  the  results  for 

the  four 

cases . 

EXTENT 

RELATIVE 

OF  IMAGE 

MAX 

MAIN  LOBE 

PSL 

ABOVE  -25  dB 

AMP 

GAIN  (dB) 

(dB) 

(mrad) _ 

DSA 

53.26 

-2.30 

-6.9 

28 

MSA 

61.34 

-1.08 

-8.4 

23 

with 

nonoverlapping  subarrays 

66.68 

-0.35 

-10.0 

18 

with 

overlapping  subarrays 

69.43 

0 

-14.4 

17 

TABLE  3.7.  COMPARISON  OF  DSA  AND  MSA,  WITH  AND  WITHOUT  SUBARRAY  PROCESSING. 

DATA  SET  309  (HIGH  STREET) .  TARGET  IS  A  CORNER  REFLECTOR.  THE 
AMPLITUDE  VARIANCE  OF  THE  BEAMFORMER  IS  0.21. 

Another  example  is  shown  in  Figure  3.20.  The  1-D  images  are 
of  a  target  in  the  Limerick  nuclear  power  plant .  The  target  is 
located  in  bin  32.  The  beamformer  is  located  in  bin  26.  The  ampli¬ 
tude  variance  of  its  echoes  across  the  array  is  0.20.  Table  3.8 
lists  the  performances  of  the  four  procedures. 


FIGURE  3.20  CROSS-RANGE  IMAGES  OF  BIN  32  WITH  LIMERICK  NUCLEAR  POWER  PLANT 
(FROM  [24]  AND  APPENDIX.) 

(a)  DSA  WITH  BIN  26 

(b)  MSA  WITH  BINS  17,  25  and  26 

(c)  MSA  WITH  NONOVERLAPPING  SUBARRAYS  USING  BINS  17,  25,  26  AND 
THREE  SUBARRAYS 

(d)  MSA  WITH  OVERLAPPING  SUBARRAYS,  USING  BINS  17,  25,  26  AND 
THREE  SUBARRAYS 
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MAX 
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13 
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-12.1 

13 

TABLE  3.8.  COMPARISON  OF  DSA  AND  MSA,  WITH  AND  WITHOUT  SUBARRAY  PROCESSING. 

DATA  SET  301  (LIMERICK  NUCLEAR  POWER  PLANT) .  THE  AMPLITUDE 
VARIANCE  OF  THE  BEAMFORMER  IS  0.20. 

The  results  in  these  two  experiments  are  similar.  Comparing 
the  top  rows  (DSA)  in  Tables  3 . 7  and  3 . 8  with  the  bottom  rows  (MSA 
with  overlapping  subarrays)  we  see  an  improvement  in  main  lobe 
gain  of  about  2  dB,  a  reduction  in  sidelobe  level  of  about  6  dB 
and  a  reduction  in  the  angular  size  of  the  target  of  nearly  2:1. 
3.10  OBSERVATIONS 

A  dominant  scatterer  algorithm  will  successfully  self-cali- 
brate  a  distorted  phased  array  when  a  single  suitable  reflector 
(large  RCS  and  small  physical  size)  is  embedded  in  clutter.  The 
target-to-clutter  ratio  must  be  large  enough  to  ensure  that  the 
reradiation  wavefront  is  dominated  by  the  prominent  reflector.  A 
test  is  the  normalized  amplitude  variance  of  the  echo  across  the 
array.  A  value  of  0.12  or  less  ensures  that  the  DSA  will  provide  a 
useful  image. 
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Under  the  same  condition,  the  MSA  will  provide  a  slightly 
superior  image  but  the  improvement  is  often  not  sufficient  to  war¬ 
rant  the  extra  complexity.  The  strength  of  MSA  is  that  it  extends 
the  useful  variance  range  to  0.2.  Subarray  processing  further  im¬ 
proves  its  performance. 
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4.  SPATIAL  CORRELATION  ALGORITHMS  (SCA) 

4 . 1  UNIT  LAG  ALGORITHM 

Figure  3.3b  showed  a  prominent  scatterer  imbedded  in  a  clut¬ 
tered  range  bin.  The  radiation  field  at  the  array  was  dominated  by 
the  large  scatterer.  As  a  consequence,  the  DSA  was  able  to  capture 
a  relatively  smooth,  linear  phasefront  and  self-calibration  was 
successful.  The  situation  pictured  in  Figure  3.3a  is  different. 
There  no  single  scatterer  is  dominant.  Instead,  the  field  at  the 
array  (3.22)  is  the  sum  of  randomly  weighted  plane  waves  arriving 
from  random  directions  within  the  illuminated  angular  sector.  Both 
situations  are  common.  Corners  figure  prominently  in  manmade  tar¬ 
gets.  Natural  terrain  such  as  farmland,  crowns  of  trees  and 
forests  is  less  structured.  Still  more  statistically  homogeneous 
is  sea  clutter.  The  spatial  correlation  algorithm  (SCA)  was  devel¬ 
oped  to  handle  these  cases  [17],  [18].  This  algorithm  is  also  due 

to  Attia. 

Let  s(u)  represent  the  complex  reflectivity  function  of  one 
of  the  range  bins  in  Figure  3.3a  s (u)  can  be  considered  a  random 
process.  For  statistically  homogeneous  clutter,  both  the  real  and 
imaginary  parts  of  s(u)  have  pdfs  that  are  independent  of  u.  The 
source  strength  is  the  random  process  s (u)  modulated  by  the  trans¬ 
mitted  beam  fT(u).  Let  S  =E{|s(u)|2}  be  the  mean  squared  value  of 

the  random  process  s  (u)  .  The  mean  square  scattered  power  as  a 
function  of  u  is  then  S|fT(u) |2. 


The  complex  scatter  strengths  in  the  range  bin  shown  in 
Figure  3.3a  represent  a  sample  function  of  the  random  process 
s(u).  It  is  reasonable  to  assume  that  the  scattering  functions 
s(u)  in  many  other  range  bins  are  also  derived  from  the  same  ran¬ 
dom  process,  especially  for  target  scenarios  such  as  sea  clutter. 
Consider  the  correlation  between  the  data  measured  at  two  array 
points  xi  and  X2,  R(xi,X2)  =  E  {e  (xi)  e*  (X2)  } ,  where  e(x)  is  the  data 
sample  measured  at  aperture  location  x.  The  Van-Cittert  Zernicke 
theorem  states  that  R(xi,X2)  is  spatially  stationary  when  the  un¬ 
derlying  source  distribution  is  incoherent,  i.e.,  R(xi,X2)  =  R(xi~ 
X2>  .  Sea  clutter  echoes  possess  this  property  because  the  source 
distribution  s  (u)  is  incoherent.  The  averaging  required  to  esti¬ 
mate  the  correlation  R(xi,X2)  can  be  performed  over  range  because 
the  source  distribution  s(u),  in  different  range  bins,  are  differ¬ 
ent  realizations  of  the  same  random  process. 

Under  far-field  conditions,  the  autocorrelation  function  R(x) 
of  the  measured  aperture  data,  and  the  mean  angular  power  density 
function  I(u)  =  S|fT(u)  1 2,  are  Fourier  transform  pairs,  as  de¬ 
scribed  by  the  Van-Cittert  Zernicke  theorem.  The  theoretical 
value  is  known  because  the  transmitted  power  pattern  is  known. 
The  measured  correlation  function  will  differ  from  the  known  theo¬ 
retical  value  because  of  phase  errors  at  the  receiving  aperture. 
The  SCA  uses  the  differences  between  the  known  and  the  measured 
correlation  phases  to  estimate  the  unknown  phase  errors  at  the  ar¬ 
ray  elements. 
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The  concept  is  introduced  in  simplified  form  by  the  following 
example.  Assume  that  |  f T  ( u )  I  is  an  even  function  of  u,  which  is 
typical  of  most  antenna  patterns.  Then  R(x)  is  real,  by  the  prop¬ 
erty  of  the  Fourier  transform.  Further,  assume  that  far-field 
conditions  prevail  and  that  the  receiving  array  is  periodic. 

These  three  assumptions,  viz.,  an  even  transmitter  pattern, 
far-field  geometries  and  periodic  receiving  array  are  not  neces¬ 
sary  for  the  SCA;  they  are  used  for  simplifying  this  disacussion. 
The  implications  of  the  the  three  assumptions  are  discussed  later. 

Figure  4.1  sketches  a  periodic  array  of  N  elements  steered  to 
0Q,  the  angle  from  broadside.  The  phase  shift  <|>n  in  the  nth  channel 

is  nkduQ  where  k  =  27C/X  is  the  wavenumber,  d  is  the 


/ 
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element  spacing,  and  uQ  =  sin0o.  There  is  a  phase  error  Pn  in  the 

nth  channel,  resulting  both  from  element  position  error  and  elec¬ 
trical  mistuning.  The  left  channel  is  considered  the  reference. 
Its  phase  error  P  is  arbitrarily  assigned  the  value  of  zero.  In 

each  channel  but  the  reference  channel  is  a  phase-error  correcting 
weight.  The  sum  is  the  array  output. 

Information  for  determining  the  correct  values  of  the  weights 
is  obtained  from  the  products  of  signals  in  adjacent  channels.  The 

It  II 

first  product  is  e2  ei*,  the  second  product  is  ^2e2*'  anc*  so  on ' 


Each  product  is  integrated  or  averaged  over  M  range  bins . 


A  T  M  I  I  £ 

R(2,1)  =  M  ^  e2,mel  'm 
m 


(4.1) 


is  an  estimate  of  the  correlation  coefficient  of  the  radiation 

field  between  the  first  and  second  antenna  elements  modified  by 
the  phase-error  difference  (32“Pi/  which,  because  Pj  is  assigned  the 

A 

value  of  zero,  equals  P2.  R(2,l)  is  a  random  variable  having  an 

expected  value 


E{R (2, 1) }  *  R (2, 1) exp ( jP2) 

A 

where  R(2,l)  is  the  error-free  correlation 


mate  of  the  correlation  coefficient  R(2,l) 
to  its  mean,  that  is 


(4.2) 

coefficient.  An  esti- 
is  approximately  equal 


R (2, 1)  ~  R(2,l)expj(p2-Pi>  (4.3a) 

In  a  similar  fashion 
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(4.3b) 


R(3,2)  ~  R (3,  2 )  exp j  ( ^3-^2 > 

R  (N,  N-l )  ~  R(N,N-l)expj  (|3n-j3n+l)  (4.3c) 

Because  of  the  far-field  assumption,  the  correlation  function 
is  spatially  stationary.  That  is,  R(i,j)  depends  only  on  (xi-xj)  . 
Further,  because  the  array  is  periodic 

R  (2, 1)  =  R (3, 2)  =  ...  =  R (N, N-l )  (4.4) 

Call  this  correlation  coefficient  R(l).  R(n+l,n)  is  real  and  posi¬ 
tive  when  (a)  both  the  transmitter  beam  and  the  scatterer  distri¬ 
bution  are  symmetrical  about  the  optical  axis  of  the  transmitting 
antenna,  and  (b)  the  spacings  between  all  adjacent  elements  are 
smaller  than  the  lobe  width  of  the  correlation  function.  In  this 
case,  the  measured  arguments  of  the  correlation  coefficients  pro¬ 
vide  estimates  of  the  elemental  phase  errors: 

p!  =0,  P2  -  Phase{R(2, 1) },  P3  =  p2  +  Phase{R(3,2) }, . .  (4.5a) 

or 

Pn  =  S”_^Phase{R(m+l,m)  },  n>l.  Pi  =  0  (4.5b) 


The  phase-compensation  weight  vector  for  the  array  is  the  conju¬ 
gate  of  the  phase-error  vector. 

4.12  UNIT  LAG  SCA  FOR  NEAR-FIELD  GEOMETRY  AND  OFF-BROADSIDE 
STEERING 

Three  constraints  were  implied  in  the  discussion  above.  First 
the  transmitter  pattern  fT(u)  was  implicitly  oriented  broadside  to 

the  receiving  array.  Second,  the  description  of  the  algorithm  was 
given  in  far-field  terms.  Third,  the  array  was  assumed  to  be  pe- 


101 


riodic.  The  implications  of  the  these  assumptions  on  the  ability 
of  the  SCA  to  calibrate  the  array  is  now  discussed. 

It  is  shown  below  that  the  near-field  assumption  is  not  nec¬ 
essary  because  the  SCA  automatically  introduces  the  phase  shift 

a 

for  focusing  the  array  into  the  phase  corrections  Pn. 

with  near-field  curvature,  the  measured  correlations,  denoted 

a 

Ri(m,n),  can  be  expressed  as 

A  A  2  2 

Rl (m, n)  ~  R(m,n)  exp [- jk -  xn) /2r ]  (4.6) 

where  xm  is  the  location  of  the  mth  element,  r  represents  range, 

/v 

R(m,n)  is  the  far-field  correlation  coefficient  given  by  (4.1)  and 
the  Fresnel  kernal  exp (- jkx2/2r)  arises  due  to  near-field  cur¬ 
vature.  The  approximation  indicates  that  the  relationship  (4.6) 
holds  in  the  region  of  validity  of  the  Fresnael  transform. 

Next,  the  effect  of  steering  away  from  broadside  is  dis¬ 
cussed.  Assume  that  uo  is  the  steering  angle.  The  function  S|fT(u 

-  uo) I2  is  the  scattered  angular  power  density  spectrum  with  the 
axis  of  symmetry  generalized  from  u  =  0  to  u  =  uo .  From  Fourier 
theory,  it  is  easy  to  show  that  the  measured  phases  of  the  corre- 

A 

lation  coefficients  R(m,n)  are  multiplied  by  exp [ jkuo (xm-xn)  ]  to 
account  for  off-broadside  steering.  With  near-field  geometry  and 
off-broadside  steering  included,  the  measured  correlation  co- 

A 

efficients,  denoted  R2(m,n),  can  be  expressed  as 

a  a  2  2  2 

R2(m,n)  -  R(m,n)  exp{  jk  [  (xm-xn)uo  -  (x^x^  (l-u0/2r]  }  (4.7) 

where  uq  is  the  steering  angle. 


102 


The  spatial  correlation  algorithm  uses  the  correlation  phases 
to  determine  the  phase  errors.  Spatial  stationarity  of  the  corre¬ 
lation  function  assumes  that  the  phase  of  the  far-field  function 

A 

R(m, n)  equals  [Pm  Pn  +  V(xm“xn)  ]  /  where  \y(xm-xn)  is  the  phase  error 
of  the  error-free  correlation  coefficient.  Therefore,  the  phase 

A 

of  R2(m,n)  in  (4.7)  can  be  expressed  as 

A  2  2 
Phase {R2 (m,n) }  =  (Pm  +  k[xmuo  -  xm(l-uQ/2r] ) 

-  (Pn  +  k[xnuo  -  x^(l-UQ/2r]} 

+  V(xm-xn)  +  8  (m,  n)  (4.8) 

The  first  (second)  term  on  the  right  is  the  sum  of  the  phase  error 
at  the  mth  (nth)  element  and  the  phase  shift  necessary  to  steer 
the  mth  (nth)  element  toward  direction  uo  and  range  r.  The  third 
term  is  the  phase  of  the  error-free  correlation  coefficient 

A 

R(m,n)  and  depends  only  on  (xm-xn)  .  The  fourth  term  8(m,n) 
represents  noise. 

The  unit  lag  SCA  uses  the  phases  of  the  (N-l)  correlation  co- 

A 

efficients  R2 (m, n)  to  determine  the  phase  errors  P2,  ...,Pn  relative 
to  Pi.  Using  (4.5b)  and  (4.8),  the  unit  lag  SCA  phase  error  esti¬ 
mates  can  be  expressed  as 

Pn  =  2-m=1Phase{R2  (m+l,m)  },  n  >1 

=  (pn  +  k[xnuo  -  xj;  (l-Uo/2r]  }  +  ^U^Cxm+l  “  xm) 


103 


^m=l  8<m+l,m) 


(4.9' 


+ 


The  first  term  on  the  right  is  the  sum  of  the  phase  error  at 
the  nth  element  and  the  phase  shift  necessary  to  steer  the  nth 
element  toward  direction  uo  and  range  r.  The  third  term  represents 
noise.  The  second  term  is  the  sum  of  the  phases  of  the  error-free 

A 

correlation  coefficient  R(m+l,m),  1  <  m  <  n-1 .  The  phase  error 

A 

estimates  (3n  become  meaningless  if  this  term  varies  randomly  as  a 
function  of  n. 

There  are  two  conditions  under  which  the  latter  term  does  not 
affect  the  performance  of  the  SCA: 

(1)  The  array  may  be  either  periodic  or  aperiodic  and  the  error- 
free  correlation  phases  ijf(m+l,m)  =  0  for  all  m,  i.e.,  the  error- 

free  unit  lag  correlation  coefficients  R(m+l,m)  are  real  and  posi¬ 
tive  . 

(2)  The  array  is  periodic. 

R(m+l,m)  will  be  real  when  the  product  of  the  transmitted 
beam  pattern  and  the  scatterer  distribution  is  symmetrical  about 
the  optical  axis.  Further,  R(m+l,m)  will  be  positive  when  all  the 
spacings  between  adjacent  elements  |  xm.1  -  xm  |  are  smaller  than 

the  null  to  null  width  of  the  correlation  function  R(x)  .  The 
width  of  the  main  lobe  of  R(x)  approximately  equals  the  length  of 
the  transmitting  antenna  Lt .  Thus  the  error-free  correlation 
phases  \|/(m+l,m)  =  0  when  (a)  the  product  of  the  transmitting  beam 
pattern  and  the  scatterer  distribution  is  symmetrical  about  the 
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optical  axis  and  (b)  |  xm_i  -  xm  |  <  Lt,  all  m.  Under  these  condi¬ 

tions,  the  phase  corrections  can  be  rewritten  as 

a  22  _ni 

Pn  =  tPn  +  kfxnuo  -  x^(l-u0/2r]}  +  Zm=1  8(m+l,m),  n>l  (4.10) 

A 

Note  that  the  phase  correction  (Sn  is  the  sum  of  the  phase  error  Pn, 
the  phase  required  to  steer  the  nth  element  to  direction  uo  and 
range  r,  and  noise. 

No  assumption  about  array  periodicity  were  made  while  deriv¬ 
ing  (4.10)  .  Thus  the  unit-  lag  SCA  is  capable  of  calibrating  ei¬ 
ther  periodic  or  aperiodic  arrays  when  (a)  the  product  of  the 
transmitter  beam  pattern  and  the  angular  power  density  spectrum  of 
the  scatterer  distribution  is  even  with  respect  to  the  optical 
axis  and  (b)  |  xm-i  -  xm |  is  smaller  than  half  the  null  to  null 

width  of  R (x) . 

4.1.3  UNIT  LAG  SCA  FOR  COMPLEX  CORRELATION  FUNCTIONS 

When  the  correlation  function  is  not  real  and  its  phase 
varies  rapidly  in  a  nonlinear  fashion  with  interelement  spacing, 
there  is  no  hope  of  using  the  correlation  measurements  to  phase  an 

aperiodic  array.  This  can  be  seen  by  observing  that  the  second 
term  in  (4.9),  <‘^m=l'^Xm+1  ~  xm) ,  will  vary  randomly  with  n.  Hence 

different  random  numbers  become  added  to  each  of  the  phase  cor¬ 
rections  Pn,  thereby  rendering  them  useless  for  phase  correcting 

the  array. 

However,  a  periodic  array  can  be  calibrated  even  under  these 
circumstances.  For  a  periodic  array,  the  second  term  in  (4.9), 
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^m=l^Xm+1  “  xm)  =  where  d  is  the  interelement  spacing 

and  \y(d)  is  the  phase  of  the  correlation  function  R(x)  for  x  =  d. 
Given  this  condition,  the  phase  corrections  in  (4.9)  can  be  ex¬ 
pressed  as 

Pn  =  <Pn  +  k  [xnuo  -  x^(l-UQ/2r]}  +  (n-l)\|/(d) 

+  8(m+l,m),  n>l  (4.11) 

In  addition  to  the  phase  errors,  the  steering  phase  and  noise,  the 
phase  correction  contains  a  linear  phase  term  (n-l)y(d).  This 

linear  phase  term  has  the  effect  of  pointing  the  array  away  from 
uo-  The  beam  pointing  error  Au  can  be  calculated  b_-  equating  the 
linear  phase  variation  (n-l)y(d)  to  k(n-l)dAu,  giving  the  beam 
pointing  error  Au  ■  \j/( d)/kd. 

Thus  periodic  arrays  can  be  phase  cohered  by  the  unit  lag  SCA 
even  when  the  phase  of  R(x)  varies  rapidly,  in  a  nonlinear  fash¬ 
ion,  with  interelement  spacing.  Such  a  phase  variation  does  cause 
a  beampointing  error.  Beampointing  errors  are  acceptable  in  ap¬ 
plications  such  as  imaging  because  the  pointing  error  leads  only 
to  an  image  shift  and  not  to  a  degradation  in  image  quality.  It 
is  valuable  also  to  note  that  when  the  phase  of  R(x)  can  be  ap¬ 
proximated  as  a  linear  function  having  slope  a  for  small  lags,  an 

aperiodic  array  can  also  be  cohered;  the  beampointing  error  in 
this  case  is  a/k. 

9.1.4  EFFECT  OF  NOISE  ON  THE  SCA 
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Eqs.(4.10)  and  (4.11)  indicate  that  the  phase  corrections 
steer  the  array  toward  uo  (or  uo  +  Au)  and  r,  in  addition  to  esti¬ 
mating  the  phase  errors.  From  (4.9),  the  noise  term  in  the  phase 
correction  for  the  nth  element,  denoted  &(>n,  can  be  expressed  as 


&|>n 


8 (m+l,m) , 


n>l 


(4.12) 


Clearly, the  number  of  noise  terms  increases  linearly  with  element 
number  n.  Thus  one  would  expect  the  variance  of  the  residual  er¬ 
rors  in  the  phase  corrections  to  increase  linearly  with  element 
number.  This  result  was  obtained  in  [17]. 

The  increase  in  the  residual  phase  error  variance  with  ele¬ 
ment  number  results  in  a  decreasing  relative  improvement  in  algo¬ 
rithm  performance  as  the  array  grows  in  size.  The  gain  of  an  ar¬ 
ray  is  expected  to  increase  directly  with  the  number  of  elements 
N.  The  rate  of  increase  with  the  SCA  is  smaller,  however. 
Because  the  phase  error  variance  grows  linearly  with  element  num¬ 
ber,  the  average  phase  error  variance  across  the  array  grows  with 
the  number  of  elements.  As  a  consequence,  the  loss  in  main  lobe 
gain  due  to  the  residual  phase  error  also  grows  with  N.  The  for¬ 
mula  for  the  normalized  main  beam  expected  gain  of  a  periodic  ar¬ 
ray  cohered  by  the  SCA  is  given  approximately  in  [17]  as 


_  .  G  _  1  2  r1/2  r  (N-l)  -  Nr  +  rN 

E[Gg]  N  N2  (1  -  r)2 


(4.13) 


where  r 


exp (-05) 


2 

and  Og  is  the  variance  of  the  noise  terms 


8 (m+1, m)  in  (4.8). 
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4 . 2  UNIT  LAG  SCA  EXPERIMENTS 


4.2.1  COMPUTER  SIMULATIONS 

In  [17],  a  CLUTTER  patch  was  simulated  on  the  computer  at  X- 
band  (X  =  3  cm)  by  uniformly  distributing  a  large  number  of 
scatterers  (1000)  within  its  boundary  (AR  =  15m,  A8  =  0.06  rad). 

The  phase  of  each  scatterer  was  taken  as  an  independent  random 
variable  uniformly  distributed  over  the  interval  [-7t,7C],  while  its 

amplitude  was  taken  to  be  proportional  to  both  the  transmitter 
radiation  pattern  and  a  uniformly  distributed  random  quantity.  The 
SCA  was  tested  through  the  simulation  of  periodic  and  random 
arrays.  Random  unknown  phase  errors  uniformly  distributed  over  the 
interval  [-71, Jt]  were  inserted  in  each  receiving  channel.  Also,  an 

independent  additive  narrowband  Gaussian  noise  component  was 
injected  into  each  channel  on  reception.  It  is  reported  in  [17] 
that  the  algorithm  was  highly  successful,  always  enabling  the 
recovery  of  the  error-free  radiation  pattern  almost  completely. 
Figure  4.3  shows  the  results  of  one  such  simulation  experiment.  It 
corresponds  to  an  N  *  20  element  random  array  with  a  nominal 
length  L  =  lOOA..  The  clutter-to-noise  ratio  (CNR)  is  20  dB  and  the 

sample  size  M  *  100.  The  solid  line  plots  the  error-free  radiation 
pattern  while  the  dotted  line  represents  the  distorted  pattern 
(due  to  the  phase  errors)  .  The  dashed  line  plots  the  pattern 
recovered  by  applying  the  SCA.  The  agreement  is  excellent. 

4.2.2  A  SEA  CLUTTER  EXPERIMENT 

In  the  Naval  Research  Laboratory  experiment  referenced  earli¬ 
er  in  Section  3.5  the  radar  data-taking  equipment  was  also  flown 
over  the  North  Atlantic  Ocean.  Eight  channels  of  sea  clutter 
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3  COMPUTER  SIMULATION  EXPERIMENT  OF  THE  UNIT-LAG  SPATIAL  CORRELA¬ 
TION  ALGORITHM.  ( - )  ERROR-FREE  PATTERN,  (...)  DISTORTED  DUE 

TO  PHASE  ERRORS,  ( - )  RECOVERED  BY  THE  ALGORITHM  (N-20, 

L-100X,  CNR- 20  dB),  M-100) .  (FROM  [21].) 
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echoes  were  recorded,  one  for  each  dipole  antenna  in  the  array,  in 
response  to  each  pulse  transmission.  The  receiving  antenna  was  a 
rigid,  8-element  equal-gain  phased  receiving  array  operating  at 
850  MHz.  Its  radiation  pattern  was  of  the  form  sin  Nx/Nsinx  where 
N  =  8  was  the  number  of  elements  and  x  is  proportional  to  the  sine 
at  the  angle  from  boresite.  Its  peak  sidelobe  is  approximately  -12 
dB.  The  pattern  is  shown  as  the  solid  curve  in  Figure  4.4. 

To  conduct  the  experiment,  random  phases  representing  ele¬ 
ment  position  errors  and  electrical  system  errors  were  added  in 
each  channel.  The  distribution  of  these  phase  errors  was  uniform 
in  the  interval  [-7C,Jt]  .  Phase  errors  destroyed  the  radiation  pat¬ 
tern,  as  expected;  it  is  shown  as  the  dotted  curve.  Sea  clutter 
echoes  measured  in  each  of  the  eight  receiving  channels  were  then 
phase-shifted  by  these  random  phases  and  applied  to  the  SCA.  100 
range  bins  of  data  were  used  to  smooth  the  correlation  estimates. 
Phase  corrections  were  calculated  from  the  measured  correlation 
coefficients  according,  to  (4.7)  and  cascaded  with  the  random  phase 
errors.  The  radiation  pattern  based  upon  the  composite  weight  vec¬ 
tor  is  shown  dashed  in  the  figure.  It  is  nearly  indistinguishable 
from  the  original,  error-free  pattern  down  to  the  third  sidelobe, 
where  the  difference  is  1.8  dB  at  the  level  of  -16  dB.  There  would 
be  no  way  to  distinguish  these  patterns  in  a  realistic  operating 
system.  This  experiment  demonstrates  the  ability  of  the  SCA  to 
phase-compensate  on  sea  clutter;  it  is  typical  of  many  similar 
tests . 


Two-dimensional  SCA  experiments  on  ground  clutter  and  on  com¬ 
plicated  echoes  from  an  airplane  also  were  performed.  These 


results  are  shown  in  Section  4.4  where  the  performances  of  several 
algorithms  are  compared. 

4.3  MULTIPLE  LAG  SCA 
4.3.1  MULLER ' S  THEOREM 

The  SCA  described  in  Section  4.1  employed  only  the  unit  lag 
correlation  coefficient  of  the  correlation  function.  The  phase  er¬ 
rors  were  derived  directly  from  the  arguments  of  the  measured  val¬ 
ues  of  R(l)  at  each  point  in  the  array.  By  the  same  procedure  that 
led  to  (4.5),  it  is  easy  to  see  that  other  correlation  coeffi¬ 
cients  contain  phase  error  information  as  well.  Consider  R(2),  for 
example.  It  provides  the  phase  error  between  two  elements,  spaced 
by  twice  the  interelement  spacing,  at  each  point  in  the  array  at 
which  it  is  measured.  Similarly,  R(3)  provides  the  phase  error  be¬ 
tween  two  elements  spaced  by  three  times  the  interelement  spacing. 
Thus  measurement  of  all  the  available  coefficients  R(l)  through 
R(N-l)  provides  phase  error  information  as  directly  as  R(l)  alone. 
The  obvious  question  to  be  answered  is  "If  R(l),  measured  every¬ 
where  in  the  array,  solves  the  phase-error  compensation  problem, 
why  utilize  any  other,  more  complicated  correlation  algorithm?" 

The  answer  is  that  three  factors  limit  the  quality  of  the 

a  A  A 

unit  lag  correlation  measurements  Rn,n+l  =  R(n,n+1)  .  The  first  is 

noise.  The  quality  of  the  estimate  is  a  function  of  the  signal-to- 
noise  ratio  (SNR)  of  the  measured  samples  of  the  electric  field,  en 

and  en+l.  Second,  the  quality  is  not  only  impaired  when  the  SNR  is 

low.  It  is  also  impaired  when  the  e-field  is  low  at  some  point  or 
points  in  the  array.  This  condition  easily  occurs  when  the  source 
function  p(u)  contains  a  small  number  of  strong  reflectors. 
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Third,  the  number  of  range  bins  used  for  smoothing  the  correlation 
coefficients  may  be  insufficient  to  produce  suitably  accurate  es¬ 
timates  of  the  phase  errors. 

When  any  or  all  of  these  conditions  prevail,  the  estimated 

A 

phase  errors,  deduced  from  the  arguments  of  R(n,n+1),  can  be  in 
error.  The  solution  is  to  measure  the  same  phase  errors  through 
other  correlation  coefficients.  In  essence,  an  additional  domain 
of  averaging  or  smoothing  is  introduced. 

Historically,  the  concept  was  first  introduced  by  Muller  as 
an  important  theorem  in  optical  telescopy  [20]  .  The  problem  of 
poor  seeing  conditions  has  plagued  astronomers  since  the  invention 
of  the  telescope.  The  problem  is  due  to  the  nonhomogeneity  of  the 
refractive  index  of  the  atmostphere.  Muller  discovered  that  it  is 
possible  to  fully  phase-compensate  a  multiple  mirror  telescope 
(analogous  to  a  phased  array  with  very  large  subapertures  for  ele¬ 
ments)  by  feedback  controlling  optical  time  delays  placed  in  front 
of  each  optical  aperture  so  as  to  maximize  the  integral  of  certain 
functions  of  the  square  of  the  image  intensity  in  the  composite 
image  plane.  Because  of  the  Fourier  relations  between  the  e-field 
in  the  aperture  and  the  image,  and  between  an  intensity  image  and 
the  autocorrelation  in  the  aperture,  Muller's  iterative  algorithm 
can  be  replaced  by  a  noniterative  function  of  the  autocorrelation 
coefficients  measured  in  the  aperture.  This  procedure  is  now  de¬ 
scribed. 

Of  the  class  of  solutions  discovered  by  Muller  the  simplest 
one  to  describe  is  as  follows.  Let  the  image  formed  by  the  ith 
mirror  be  1^  Let  all  the  images  be  brought  together  onto  a  common 
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focal  plane  and  registered.  Form  an  objective  function  J,  called 
the  sharpness  function,  that  is  the  integral  of  the  square  of  the 
intensity  distribution  in  the  image  plane.  Let  the  time  delays  re¬ 
ferred  to  above  be  represented  more  generally  by  a  set  of  weights 
Wj_.  The  optimization  or  search  procedure  adjusts  the  w^  until 

9j 

- —  =  0,  all  i  (4 . 14) 

OW  2. 

The  theorem  demonstrates  that  when  this  condition  is  achieved  the 
wavefront  distortion  has  been  fully  compensated  and  that  the  image 
is  without  error,  except  for  an  arbitrary  and  unknown  displacement 
of  the  optical  axis. 

In  mathematical  terms,  let  s  =  s  (u)  be  the  source  distribu¬ 
tion,  f Q~f Q ( u )  the  design  radiation  pattern  and  f  »  f(u)  «  fQ  +  8f 

the  actual  pattern.  The  optimum  complex  image  is 

sQ  =  s  *  fQ  (4.15) 

while  the  actual  distorted  image  is 

s  =  s*f  =  s*  (fQ  +  8f) 

-  sc  +  8s  (4.16) 

The  display  I  =  I(u)  is  either  the  magnitude  or  the  square  magni¬ 
tude  of  (4.16)  and  IQ  is  displayed  when  the  radiation  pattern  is 

error  free .  The  statement  of  the  theorem  is  that 

I(u)  =  IQ(u  +  8u)  (4.17) 
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when  (4.14)  is  satisfied,  where  8u  is  the  shift  in  the  optical 
axis . 

A  limitation  of  the  theorem  is  the  requirement  that  the  radi¬ 
ation  field  upon  which  it  operates  be  derived  from  an  incoherent 
distribution  of  sources.  Stars,  being  independent  energy  genera¬ 
tors,  satisfy  this  requirement.  The  phase  relations  between  the 
radiations  from  each  source  within  the  passbands  of  the  receivers 
are  continuously  changing.  The  result  is  a  radiation  field  at  the 
receiving  antenna  which  is  continuously  changing  but  whose  statis¬ 
tics  are  invariant  with  time. 

4.3.2  THE  MICROWAVE  CASE 

Not  so  in  microwave  radar,  however.  There  the  phase  relations 
between  the  echoes  from  targets  are  dictated  by  the  geometry  of 
the  locations  of  the  transmitter,  the  targets,  and  the  receiver. 
Each  transmitted  pulse  results  in  the  same  set  of  echoes  and  hence 
the  same  radiation  field  at  the  receiving  array. 

Figure  4.5  schematically  portrays  a  surface  J  as  a  function 
of  two  weights,  indicated  by  the  phase  shifts  <|)1  and  <J>2  when  the 

radiation  field  is  not  incoherent.  The  surface  shown  has  multiple 
maxima  and  the  search  procedure  can  mistakenly  settle  for  a  non- 
global  peak.  The  optimum  weights  are  the  coordinates  of  the  maxi¬ 
mum  modal  point  of  the  surface . 

The  primary  potential  benefit  of  the  procedure,  in  microwave 
radar,  is  less  its  ability  to  compensate  for  medium  turbulence 
than  for  distortion  in  the  aperture.  Large  antenna  systems  will  be 
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intrinsically  distorted  due  to  nonrigidity  of  their  structures, 
differences  in  gravitational  forces  across  the  antennas,  differen¬ 
tial  windloading,  differential  weather,  etc.  To  achieve  the  bene¬ 
fits  of  self-calibration  based  upon  Muller's  theorem,  a  means  must 
be  developed  for  altering  the  phase  relations  of  the  electric 
fields  reflected  from  the  targets  so  as  to  emulate  an  incoherent 
distribution  of  radio  sources. 

Two  methods  for  doing  so  have  been  developed.  Each  method 
builds  up  a  diversity  image  suitable  for  initiation  of  the  search 
for  the  weight  vector.  The  first  procedure  uses  transmitter  loca¬ 
tion  diversity  as  its  tool  [19].  The  second  procedure  uses  range 
bin  diversity  (Section  4.1).  Transmitter  diversity  solves  the 
problem  completely.  Its  limitation  is  its  practicality,  however, 

for  it  is  not  a  solution  applicable  to  all  system  designs.  The  di¬ 
versity  image  ID  is  the  weighted  sum  of  images  Im,  m=l, 2,  . . .M,  ob¬ 
tained  from  a  set  of  M  snapshots  of  the  radiation  field.  These  im¬ 
ages  are  called  trial,  images  or  constituent  images.  The  trial  im¬ 
ages  differ  either  because  the  diversity  procedure  alters  the  ra¬ 
diation  pattern  between  images  or  because  the  procedure  alters  the 
phase  relations  among  scatterers  in  the  source  function,  or  both. 
Thus  the  diversity  image  to  which  the  search  procedure  (3.57)  is 
applied  can  be  written 

Id  “  £  anIB  (4-18) 

m=l 

9.3.3  RELATIONSHIP  BETWEEN  MULLER'S  OBJECTIVE  FUNCTION  AND 
APERTURE  CORRELATION  FUNCTIONS 
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Assuming  that  the  characteristics  of  the  diversity  image 
(4.18)  are  consistent  with  an  incoherent  distribution  of  sources, 
its  inverse  Fourier  transform  is  the  autocorrelation  function  R  = 
R ( x )  in  the  aperture,  i.e., 

R  =  F-1I  (4.19) 

and  the  inverse  Fourier  transform  of  I2  is 

F-1I2  =  F-1I  *  F-1I  =  R  *  R  (4.20) 

Muller's  objective  function 

J  =  Jl2 (u) du  =  R (x) *R (x)  |  X=Q  =  Jr2 (x) dx  (4.21) 

or  in  discrete  terms 

N-l 

J=  |R(0)  I2  +  2  ^\R(1)  I2  (4.22a) 

1=1 

where 

N-i 

R(2)  =  Ajexp  ( jxj/j  )  =  £R(n+2,n),  i  >  0  (4.22b) 

n=l 

R(-i)  =  R*  (2)  (4.22 

Thus  Muller's  objective  function  can  be  calculated  from  cor¬ 
relation  coefficients  measured  in  the  aperture.  An  iterative  pro¬ 
cedures  developed  by  Subbaram  is  now  described.  This  algorithm 
uses  the  complete  correlation  function  rather  that  its  phase 
alone.  The  advantage  of  using  the  complete  correlation  coeffi¬ 
cient  lies  in  the  accuracy  of  the  phase  measurement.  When  the  mod¬ 
ulus  is  large  the  phase  is  easy  to  measure,  in  which  case  both 
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procedures  are  equally  effective.  When  the  modulus  is  small,  how¬ 
ever,  the  phase  measurement  is  subject  to  considerable  error  due 
to  noise  and  the  algorithm  suffers.  By  using  the  entire  complex 

A 

coefficient  the  amplitude  of  R  acts  as  a  weight  upon  the  measured 
phase.  In  this  way  noise-induced  phase  measurement  errors  are  sup¬ 
pressed. 

The  measured  values  of  R(n+l,n)  will  be  in  error  because  of 
aperture  and/or  medium  distortion.  The  corrolation  coefficient  for 

A 

lag  1,  after  introducing  the  phase  corrections  Pn^  may  be  written 


R1 


N-l  a 

I  R (n+1, n) 
n=l 


exP [ j (Pn  "  Pn+1> ] 


(4.23) 


Referring  to  (4.22),  observe  that  R(0)  is  independent  of  both 
the  phase  errors  and  phase  corrections.  Thus  it  can  be  dropped 
from  the  objective  function.  Further,  the  correlation  coeffi¬ 
cients  for  large  lags  have  a  tendancy  to  be  noisy;  therefore  it  is 
better  to  use  only  a  limited  number  of  correlation  lags  in  the  ob¬ 
jective  function.  Modifying  J  in  (4.22a)  by  using  only  lags  one 
through  L  results  in 


F 


L 

I 

1=1 


N-l  A 

X  R (n+1, n)  exp[-j(pn+2 
n=l 


Pn>J 


=  Fi  +  F2  + 


(4.24) 


The  first  term  is 


Fl 


N-l  A 

X  R (n+1, n)  exp[-j(pn+1 
n=l 


(4.25) 
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An  algorithm  that  maximizes  Fi  with  respect  to  the  phase  correc- 

A 

tions  Pn  may  be  called  a  unit  lag  spatial  correlation  algorithm. 
Maximizing  the  additional  terms  in  F  generalizes  the  spatial  cor¬ 
relation  algorithm  up  to  the  Lth  lag  by  attempting  to  cophase  the 
measured  correlation  data  from  lags  1  through  L.  When  the  correla¬ 
tion  data  are  distorted  only  by  phase  errors,  i.e., 

R  (n  +  1,  n)  =  Aj  exp[j(\j/i  +  (5n+i  -  (5n)  ] ,  (4.2  6) 

A  A 

the  phase  corrections  (3n  will  cophase  all  the  R(n+2,n)  for  lags  one 

A 

through  L.  In  other  words,  the  phase  correction  (3n  will  make  F 
equal  to  its  global  maximum  when  the  measured  correlation  data 

A 

R(n+l,n)  are  error-free  except  for  a  phase  shift  of  Pn+j  -  Pn-  In 
realistic  situations,  the  amplitudes  of  the  correlation  data  hav¬ 
ing  the  same  lag  1  will  not  be  constant  because  of  noise  and  an 
insufficient  number  of  range  bins  available  for  smoothing  the  cor¬ 
relation  data.  Further,  the  phases  of  the  measured  correlation 
data  will  not  be  equal  to  \| fi  +  (3n+j  -  (3n  due  to  the  same  reason. 

Hence  F  cannot  reach  its  global  maximum  when  the  measured  correla¬ 
tion  data  are  contaminated  by  noise.  The  set  of  phase  corrections 
obtained  by  maximizing  F  will  attempt  to  cophase  the  correlation 
data  giving  the  largest  weight  to  the  lags  that  have  the  largest 
amplitude,  and  small  weights  to  lags  that  have  small  amplitudes. 
Thus  maximizing  F  in  order  to  obtain  the  phase  corrections  over¬ 
comes  the  three  deficiencies  of  the  spatial  correlation  algorithm 
mentioned  at  the  beginning  of  this  section. 

9.3.4  MULLER ' S  THEOREM  IN  THE  APERTURE  DOMAIN 


121 


The  following  proposition  derives  an  expression  for  the  phase 

A 

corrections  Pn  that  are  obtained  by  maximizing  F.  A  computationally 
efficient  algorithm  for  maximizing  F  is  then  presented. 

PROPOSITION 

Maximizing  the  objective  function  F  with  respect  to  the  phase 

A 

corrections  Pn  yields  the  following  set  of  phase  corrections: 

Pn  =  Pn  +  cD  +  (n-Dc^  1  <  n  <  N  (4.27) 

A 

In  other  words,  the  phase  corrections  pn  obtained  by  maximiz¬ 
ing  F  are  equal  to  the  sum  of  the  phase  errors  pn,  an  arbitrary 
constant  cQ,  and  an  arbitrary  phase  tilt  having  a  slope  ci .  The 
constant  cQ  is  of  no  importance  because  it  merely  adds  a  phase  c0 
to  each  image  pixel.  The  constant  Cj  shifts  the  image  by  an  amount 
Au  =  -ci/kd  for  an  array  with  interelement  spacing  d.  Thus  the 

A 

phase-corrected  image  s(u)  is  a  version  of  the  error-free  image 
so(u)  multiplied  by  a  constant  eico. 

PROOF  OF  PROPOSITION 

The  phase  corrections  obtained  by  maximizing  F  are  shown  to 

have  the  form  given  by  (4.27).  Any  set  of  N  quantities  ao/.../3N-i 

can  be  written  in  terms  of  a  polynomial  of  order  N-l,  i.e.,  an  = 
N-l 

]^C)c  (n-l ) k,  where  the  c*  are  the  coefficients  of  the  polynomial. 
k=0 

A 

Expanding  the  Pn  -  Pn  in  this  fashion  leads  to 


Further,  assume  that  the  only  distortion  in  the  measured  cor¬ 
relation  data  is  due  to  the  phase  errors  (3n,  as  in  (4.26). 

Clearly,  F,  defined  in  (4 . 24 ), reaches  its  maximum  when  all  compo¬ 
nents  in  the  summation  over  n  are  cophased.  Consider  a  single 
component  Fj  of  F.  From  (4.24) 


F  i  = 


N-2  a 

X  R  (n+1,  n)  exp[-j(pn+I 
n=l 


(4.29) 


Substituting  (4.28)  and  (4.26)  into  (4.29)  yields 

N-l 

Fi  -  Aj  |  Xexp[-j  (Pn  +  i  -  pn+2  -  Pn  +  Pn>  3  I  2 

n=l 


N-l 

I 

n=l 


N-l 


exP  {  3  Xck  [  (n-l)  k-  (n+2-1)  k] 
k=o 


(4.30) 


Fi  reaches  its  maximum  when  each  term  in  the  outer  summation 
over  n  has  the  same  phase.  This  occurs  when  the  exponential  terms 
in  the  summation  are  independent  of  n.  Consider  a  single  exponen¬ 
tial  term  inside  the  summation.  Its  exponment  is  the  sum  of  N 
terms  cic[(n-l)k  -  (n+1-1) k],  0  <  k  <  N-l.  For  k  =  0,  this  quantity 
is  co  (1-1)  =  0,  which  is  independent  of  n.  The  second  or  k  =  1 
term  is  ci  (n-l-n-1-1)  =  cj.1,  which  is  a  function  of  1  but  not  of  n. 
The  third  term  is  cz(21  -2nl  -l2).  This  term  is  a  function  of  n. 
Therefore  the  phase  of  each  term  in  the  summation  over  n  will  be 
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Therefore  the  phase  of  each  term  in  the  summation  over  n  will  be 
different  if  02  *  0 .  The  requirement  that  each  term  in  ther  summa¬ 
tion  be  cophased  at  the  maximum  of  F i  therefore  demands  that  C2  = 
0.  In  a  similar  fashion,  it  is  easy  to  show  that  all  the  c*,  k  > 
2,  are  zero  at  the  maximum  Fj.  The  same  conclusion  holds  for  all 
components  of  F,  i.e.,  for  all  Fj,  \  £  1  £  L.  Using  c*  =  0,  k  >  2, 
in  (4.28)  yields 

Pn  =  Pn  +  Cc  +  (n-l)c1  (4.31) 


as  claimed  in  (4.27). 

The  following  discussion  shows  that  the  linear  phase  error 
(n-l)ci  shifts  the  image  from  the  origin  by  an  amount  Au  =  -ci/kd, 

where  d  is  the  interelement  spacing.  This  can  be  seen  by  examin- 
ing  the  phase-corrected  image  s(u)  given  by 


s  (u) 


A 

1  -jkdu(n-l)  e-jPn 


(4.32) 


The  en  are  equal  to  the  error-free  data  samples  en  shifted  in 

A 

phase  by  the  phase  error  pn.  Using  en  *  enexp(jPn)  in  (4.32)  leads 


A 

s(u)  =  I  en  e"j(pn  '  ^  e-^du(n-l)  (4.33) 

n-1 

Substitution  of  (4.31)  into  (4.33)  yields 

s(u)  .  e-jc°  I  en  e-^d(n_1)  (u+ci/kd) 
n-1 
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e"jc°  sQ(u+c1/kd) 


(4.34) 


Thus  the  phase-corrected  image  s(u)  =  is  a  shifted  version  of  the 
error-free  image  sQ(u)  multiplied  by  a  constant  exp(jcQ). 

9.3.5  A  COMPUTATIONALLY  EFFICIENT  ALGORITHM  FOR  MAXIMIZING  THE 
OBJECTIVE  FUNCTION  F 

The  objective  function  F  has  to  be  maximized  with  respect  to 

A  A  A  A 

(N-l)  parameters  P2,  P3  ...'  Pn  (remember  that  =  0).  Setting 

A 

the  derivatives  of  F  with  respect  to  the  (Pn)  to  zero  yields  a  set 

of  (N-l)  transcendental  equations  for  the  (N-l)  unknowns.  No 
closed  form  solution  for  this  set  of  equations  has  been  found; 
hence  F  has  to  be  maximized  iteratively. 

There  are  four  major  steps  in  the  following  iterative  algo¬ 
rithm: 

A  A 

1.  Choose  P1  *  0  and  Pn,  n  >  1,  using  the  spatial  correla¬ 
tion  algorithm  (Section  4.1). 

*  A 

2.  Maximize  F  sequentially  with  respect  to  the  P^,  2  <  k  < 

A 

N,  while  keeping  the  rest  of  the  pn,  n*  k,  constant. 
Calculate  the  global  maximum  Fmax  from  the  data. 

3.  Check  if  F  is  close  to  Fmax .  If  not,  repeat  step  2. 

Otherwise,  go  to  step  4. 

4.  Modify  the  set  of  phase  corrections  to  minimize  the  dis¬ 
placement  of  the  optical  axis. 

Step  1  is  described  in  Section  4.1.  The  remaining  steps  are 
now  described. 


125 


Step  2.  Because  F  has  several  local  maxima,  conjugate  gradi¬ 
ent  algorithms  cannot  be  used  to  accomplish  this  task.  An  algo¬ 
rithm  that  maximizes  F  efficiently  is  now  developed.  This  algo¬ 
rithm  exploits  the  fact  that  F  is  almost  a  sinusoidal  function  of 

A  A 

a  single  Pk  when  the  other  Pn,  n*k,  are  held  constant,  that  is 

A  A 

F(pk)  ~  Aksin  (p  k  +  <)>k)  +  Bk  (4.35) 

where  Ak,  <f>k  and  Bk  are  unknown  constants.  To  prove  (4.35),  con- 

A  A  A 

sider  a  single  component  of  F  (Pk)  (say  F ( Pk )  )  in  terms  of  Pk 

A 

holding  the  Pn,  n*k,  constant. 

A 

In  (4.29)  the  terms  that  contain  Pk  are 

a 2  =  R (k, k+1) exp (- jpk+1)  if  k+1  <  N 

=  0  otherwise 

and 

A  A 

a£  =  R  (k,  k+1)  exp  (  jPk_j_)  if  k-1  >  1 
=  0  otherwise 
The  remainder 

a3  =  I  R(n,n+l)exp(-j  (pn+1  -  pn)  ] 
n=l 

n  *  k  or  n  *  k-1 

A 

is  independent  of  Pk.  Thus 

Fl(Pk)  -  +  a2e"^k  +  a3  2  (4.36) 

Because  the  algorithm  begins  with  the  unit  lag  spatial  correlation 
method  (step  (1))  which  provides  high  quality  initial  estimates  of 
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the  Pn,  all  components  of  the  sum  a3  are  nearly  cophased  and,  be¬ 
cause  33  is  the  sum  of  N-l-2  terms  each  of  the  type  a^  or  83,  1 33 1 
>>  |  a  2  |  and  |  a.2 1  •  Using  this  fact  in  (4.36)  yields  an  expression 

A 

for  FifPfc)  of  the  form 

Fx(Pk)  ~  A  sin  (pk  +  <{»  +  B,  1  <  1  <  L 


where  terms  containing  exp(±j2Pk)  have  been  ignored  because  their 
magnitude  1 3^3 1  is  much  smaller  than  the  magnitude  of  the  terms 

A 

containing  exp(±jPk).  The  sum  of  1  sinusoids  of  the  form  shown 

A 

above  is  also  a  sinusoid;  hence  F  (Pk  )  is  of  the  form  shown  in 
(4.35)  . 

A. 

The  iterative  algorithm  first  maximizes  F  with  respect  to  P3, 

A 

then  does  the  same  for  p3,  etc.  After  going  through  a  cycle  of 

A 

maximizing  F  with  respect  to  the  pn,  2  <  n  <  N,  the  value  of  F  is 
compared  to  Fmax .  The  iterations  are  stopped  when  F  is  close  to 
Fmax-  Equation  (4.35)  is  used  in  order  to  maximize  F  with  respect 

A  A 

to  a  single  Pk  while  holding  the  remaining  Pn,  n*k,  constant.  Note 

A  A 

that  F(Pk)  reaches  its  maximum  value  (=  Ak+Bk)  when  Pk  =  (tc/2)  -<t>k- 

A 

The  value  of  <j>k  can  be  computed  by  evaluating  F(pk)  at  three 
points.  For  example, 


and 


=  F(0)  *  Ak  sin  ^k  +  Bk 

A  *  c 

S3  *  F (J)  =  Ak  COS  0k  +  bk 

A  .  . 

S3  -  F  (it)  =  -Ak  sin  <})k  +  Bk 
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Then  tan 


h  is 


4>k 


(S1-S3) / (2S2“Si~S3) ,  and  hence  the  optimal  value  of 


Pk 


opt 


Jt 

2 


tan-* 


t — _ 

'2S2-S1-S3 


(4.37) 


Thus  maximizing  F  with  respect  to  a  single  pk  involves  computing  F 

A 

for  only  three  values  of  (3k  . 

Step  3.  Any  iterative  algorithm  needs  a  stopping  criterion. 
Fortunately,  the  value  of  the  global  maximum  of  F  can  be  computed 
from  the  data  because  each  term  in  the  summation  over  n  in  (4.24) 

will  have  the  same  phase  when  F  reaches  its  maximum.  Thus,  the 
maximum  value  of  F,  denoted  Fmax,  i.e., 


Fmax  -  £ 

1*1 


N-I 

X 


n=l 


Rn+l,n | 


(4.38) 


Step  4 .  It  is  useful  to  minimize  the  shift  of  the  phase-cor¬ 
rected  image  (=  -c^/Rd  from  (4.34)).  Clearly,  the  magnitude  of  the 

phase-corrected  image  will  be  identical  to  the  magnitude  of  the 
error-free  image  if  c1  ■  0.  The  shift  of  the  phase-corrected  image 
will  be  minimized  if  c^  can  be  estimated. 

Consider  the  phase-corrected  unit  lag  phase;  that  is,  con¬ 
sider 


R(n+l,n) 


Vi  "  ci 


(4.39) 
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The  identity  in  (4.39)  results  when  (4.26)  and  (4.31)  are  used  for 

A  A 

R(n,n+1)  and  (3n,  respectively.  Note  that  ci  =  0  if  the  phase  of  the 


corrected  correlation  data  for  unit  lag  are  set  to  \|/i .  Hence  ci  can 


be  set  to  zero  if  the  unit  lag  correlation  phase  is  known.  For 


most  target  scenarios,  the  unit  lag  autocorrelation  phase  is  ap¬ 


proximately  equal  to  -kduQ,  where  d  is  the  interelement  spacing  and 


uQ  is  the  center  of  the  field  of  view  (the  direction  towards  which 


the  transmitter  is  pointed)  .  Therefore,  an  estimate  of  \j/i,  denoted 

A 

\\fi,  is  -kduQ.  Thus  ci  can  be  estimated  by  comparing  the  phase  of 


the  unit  lag  phase-corrected  correlation  values  to  i.e.. 


A  A  f  A  3  ($  ,  T  (3  )  •) 

ci  =  Vi  “  Phase  \  R (n+1, n)  e  n+i  n  /  (4.40) 


Thus  the  shift  in  the  image  plane  can  be  minimized  after  maximiz- 

A  A 

ing  F  by  subtracting  out  (n-l)ci  from  the  Pn .  The  final  set  of 

A  I 

phase  corrections,  denoted  pn,  are 


Pn  -  Pn  -  (n-l)C!  (4.41) 


where  the  Pn  are  obtained  by  maximizing  F  and  ci  is  defined  by 


(4.40) . 


4.4  MULTIPLE  LAG  SIMULATION  EXPERIMENTS 


Figure  4.6  is  a  computer  simulation  experiment  in  which 


transmitter  diversity  performed  the  decohering  function.  Two  equal 


strength  point  targets  are  located  at  u  =  -0.3  and  0.5.  The  solid 


curve  shows  the  one-dimensional  image  of  the  two  targets  obtained 


using  an  11-element  periodic  receiving  array.  The  dotted  curve 
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shows  the  useless  image  that  resulted  when  random  phase  errors 
were  added  at  each  receiver.  The  error  distribution  was  uniform 
over  the  interval  +  Jt. 

Images  were  then  computed  for  five  different  locations  of  the 
transmitter.  The  transmitter  locations  were  confined  to  the  extent 
of  the  receive  array.  Cross  products  obtained  from  the  five  sets 
of  receiver  measurements  were  averaged  to  provide  the  estimates  of 
the  correlation  coefficients,  which  were  then  applied  to  the  algo¬ 
rithm.  The  result  is  shown  in  the  dashed  curve  which  is  a  nearly 
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U«sin6 

FIGURE  4.6  COMPUTER  SIMULATION  EXPERIMENT  TO  TEST  MULTIPLE-LAG 
CORRELATION  ALGORITHM. 

(a)  IMAGE  WITH  NO  PHASE  ERRORS:  (b)  IMAGE  WITH  PHASE 
ERRORS:  (c)  IMAGE  WITH  PHASE  ERRORS  CORRECTED  BY  THE 
MULTIPLE-LAG  CORRELATION  ALGORITHM.  (FROM  [23].) 
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identical  copy  of  the  original  solid  curve  displaced  0.68  units  of 
u  to  the  right . 

Figure  4.7  shows  another  experiment.  In  this  experiment  a 
continuous  target  from  n1  =  -  0.55  to  u2  =  -  0.15  having  a  random 

rough  surface  was  simulated  by  one-hundred  unit  amplitude  and  ran¬ 
dom  phase  scatterers  randomly  distributed  from  ux  to  u2 .  A  combina¬ 
tion  of  range  bin  diversity  and  transmitter  diversity  was  em¬ 
ployed.  Two  range  bins  were  located  1300  X  and  1340  X  from  the  re¬ 
ceiving  array.  The  scatterer  distributions  in  the  two  range  bins 
were  independent  of  each  other.  The  receiving  array  consisted  of 
24  elements  spaced  one  half  wavelength  apart.  The  FOV  was  illumi¬ 
nated  from  three  locations,  the  transmitters  being  placed  at  the 
center  and  the  two  ends  of  the  array.  Zero  mean  Gaussian  noise 
having  a  variance  of  0.01  was  introduced  at  each  receiver.  The 
combination  of  three  transmitters  and  two  range  bins  provided  six¬ 
fold  diversity  for  decorrelating  the  target  scenario. 

Curve  (a)  in  Figure  4.7  is  the  error-free  image  of  the  first 
range  bin.  Independent  phase  errors,  uniformly  distributed  between 
-K  and  71,  were  then  introduced  at  each  receiver.  Curve  (b)  is  the 

image  of  the  first  range  bin  without  correcting  for  the  phase  er¬ 
rors.  This  image  does  not  resemble  the  target  scenario  whatsoever. 
Curve  (c)  is  the  image  of  the  first  range  bin  after  correcting  for 
the  phase  errors.  This  image  is  almost  identical  to  the  error-free 
image  except  for  a  shift  Au  =  1/3  in  the  image  plane. 

4 . 5  COMPARISONS  WITH  OTHER  ALGORITHMS 
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INTENSITY 


(a)  ERROR  FREE  IMAGE 

(b)  IMACE  WITH  PHASE  ERRORS 

(c)  IMAGE  AFTER  CALIBRATION 


CROSSRANGE  (  u  ) 


FIGURE  4.7 


DEMONSTRATION  OF  THE  MULTIPLE  LAG  CALIBRATOR  FOR  FAR-FIELD  TARGETS 
100  COHERENT  TARGETS,  LOCATED  IN  THE  FAR  FIELD  OF  THE  ARRAY  BETWEEN 

^  fnc*  u””®*15»  ARE  IMAGED  USING  a  24-ELEMENT  FILLED  PERIODIC  ARRAY. 
CURVE  (a)  IS  THE  ERROR-FREE  IMAGE  OF  THE  TARGET  SCENARIO.  CURVE  (b) 

IS  THE  IMAGE  OBTAINED  AFTER  ADDING  RANDOM  PHASE  ERRORS  AT  EACH  RECEIVER. 
CURVE  (c)  IS  THE  IMAGE  OBTAINED  AFTER  CORRECTING  FOR  THESE  PHASE  ERRORS* 
USING  THE  FIRST  FIVE  LAGS.  SIXTH  ORDER  TRANSMITTER  LOCATION  DIVERSITY  WAS 
USED  FOR  DECORRELATING  THE  COHERENT  TARGETS.  OBSERVE  THAT  THE  IMAGES  (a) 
AND  (c)  ARE  ALMOST  IDENTICAL  EXCEPT  FOR  A  SHIFT  Au-1/3  ’u'  UNITS  IN  THE 
IMAGE  PLANE. 
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The  experiments  reported  in  this  section  fall  into  two 
groups.  In  the  first  three  tests  (Figures  4.8-4.10)  the  targets 
satisfied  the  dominant  scatter  condition;  in  the  remainder  the  low 
variance  beamformer  condition  was  not  met.  The  upper  image  of  Fig¬ 
ure  4.8  was  processed  using  the  dominant  scatterer  algorithm;  tha 
lower  with  the  modified  Muller-Buff ington  procedure.  The  aircraft 
was  a  large  jet  at  a  distance  of  about  20  km  from  the  Valley  Forge 
laboratory.  Because  the  dominant  scatterer  condition  was  satisfied 
a  satisfactory  adaptively  beamformed  image  was  obtained.  The  bot¬ 
tom  figure  shows  the  image  obtained  using  the  same  data  set  when 
the  MMB  algorithm  was  applied.  It  is  evident  that  the  two  proce¬ 
dures  were  equivalent . 

Figure  4.9  compares  images  of  a  Boeing  727  formed  by  four  al¬ 
gorithms,  the  DSA,  MSA,  SCA  and  MMB.  A  reference  scatterer  with 
very  low  normalized  amplitude  variance  (O2  =  0.05)  was  on  the  tar¬ 
get.  All  methods  work  well.  The  main  lobe  gains  of  the  adaptively 
formed  arrays  differ  by  no  more  than  0.5  dB. 

Figure  4.10  compares  images  of  the  Limerick  nuclear  power 
plant  obtained  from  the  same  four  algorithms.  Again  the  results 
are  nearly  alike.  The  array  gains  differ  by  less  than  0.5  dB. 

The  situation  changes  when  no  low  variance  range  bin  can  be 
found.  Figure  4.11  shows  the  use  of  bin  47  for  which  O2  =  0.17  is 

the  reference  bin.  This  high  variance  bin  resulted  in  failure  by 
both  DSA  and  MSA;  the  former  is  shown  in  plate  a.  SCA  and  MMB  were 
both  able  to  form  images  recognizable  as  airplanes. 

Four  images  from  another  data  set  obtained  from  the  same  air¬ 
plane  are  shown  in  Figure  4.12.  Again  the  SCA  and  MMB  images  are 
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RANGE  (MEIERS)  RANGE  (METERS) 

10.0  20.0  30.0  40.0  50.0  0.0  10.0  20.0  30.0  40.0  50.0 


FIGURE  4.8  RADAR  CAMERA  ISAR  IMAGES  OF  A  LARGE  JET  AIRCRAFT  AT 

A  DISTANCE  OF  20  km  FROM  VALLEY  FORGE  RESEARCH  CENTER. 
THE  RANGE  CELL  IS  3  METERS,  (a)  DSA:  (b)  MULTIPLE  LAG 
SCA.  BOTH  PROCEDURES  WORKED  WELL  AND  YIELDED  COMPARABLE 
IMAGES.  (FROM  125].) 
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FIGURE  4.9  :  RANGE-AZIMUTH  IMAGE  OF  BOEING  727  AIRCRAFT  OBTAINED  USING  THE 
FOURTH  QUADRANT  (ELEMENTS  FROM  385  THROUGH  512).  (FROM  [26].) 

(a)  DSA  WITH  BIN  40  HAVING  AMPLITUDE  VARIANCE  0.05  AS  THE 
REFERENCE  RANGE  BIN 

(b)  MSA  WITH  BINS  29,  40  AND  42 

(c)  SCA  WITH  LAG  ONE  AND  BINS  37  THROUGH  40 

(d)  MMB  WITH  LAG  (1,6,28)  AND  BINS  37  THROUGH  40 
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FIGURE  4.11  RANGE-AZIMUTH  IMAGE  OF  BOEING  727  AIRCRAFT.  (FROM  [26].) 

(a)  DSA  WITH  BIN  47  HAVING  AMPLITUDE  VARIANCE  0.17 
AS  THE  REFERENCE  RANGE  BIN 

(b)  MSA  IS  NOT  AVAILABLE  WITH  BINS  44  THROUGH  47 

(c)  SCA  WITH  LAG  ONE  AND  BINS  44  THROUGH  47 

(d)  MMB  WITH  LAG  (1,6,28)  AND  BINS  44  THROUGH  47 
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FIGURE  4.12:  RANGE- AZIMUTH  IMAGE  OF  BOEING  727  AIRCRAFT  OBTAINED  USING  THE 
SECOND  QUADRANT  (ELEMENTS  FROM  129  THROUGH  256)  .  (FROM  [26].) 

(a)  DSA  WITH  BIN  40  HAVING  AMPLITUDE  VARIANCE  0.15  AS  THE 
REFERENCE  RANGE  BIN 

(b)  MSA  WITH  BINS  38.  39  AND  40 

(c)  SCA  WITH  LAG  ONE  AND  BINS  39  THROUGH  42 

(d)  MMB  WITH  LAG  (1,6,28)  AND  BINS  39  THROUGH  42 


superior.  Main  lobe  gain  is  3.3  dB  greater  than  DSA  and  3.9  dB 
greater  than  MSA. 

Figure  4.13  shows  DSA  and  SCA  1-D  images  of  a  corner  reflec¬ 
tor  in  which  echoes  from  farmland  were  used  for  self-calibration. 
DSA  fails  completely;  the  DSA  image  bears  no  relation  to  that  of  a 
corner  reflector  while  the  SCA  image  shows  a  prominent  target.  MMB 
performance  is  comparable  to  SCA.  Both  exhibit  4.5  dB  greater  main 
lobe  gain. 

The  observations  from  this  and  the  preceding  chapter  lead  to 
the  following  conclusions: 

•  If  at  least  one  range  bin  has  a  low  variance  beamformer  (C2 
<  0.12),  any  algorithm  is  satisfactory.  In  general  the  DSA 
is  the  simplest  to  implement. 

•  If  no  low  variance  bin  can  be  found  but  one  or  more  bins 
exhibit  single  targets  with  lower  ratios  of  target  to  clut¬ 
ter  strengths,  then 

DSA  will  fail. 

MSA  will  often  succeed,  particularly  with  the  use  of 

overlapping  subarrays . 

SCA  and  MMB  will  always  succeed. 

•  If  no  bins  exhibit  single  targets  in  clutter 

DSA  and  MSA  will  fail. 

SCA  and  MMB  will  often  succeed. 
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5.  RELATIVE  COMPUTATION  TIMES 

Comparisons  of  CPU  times  were  made  for  the  several  algorithms 
for  two  target  scenes,  the  Limerick  Nuclear  Power  Plant  and  a 
Boeing  727  ISAR  image.  The  main  distinction  between  them  is  that 
the  former  scene  is  echo-filled  in  range  and  azimuth  whereas  the 
latter  target  is  in  the  clear  and  the  scene  is  devoid  of  clutter. 
Taking  the  CPU  time  required  for  DSA  to  produce  a  1-D  image  as 
unity,  a  summary  of  the  relative  times  is  given  in  Table  5.1.  The 
reader  is  referred  to  Section  5.2  of  the  Appendix  for  further 
details . 

_ IdQ _  _ z=n _ 


a 

b 

c 

a 

b 

c 

DSA 

1 

1 

28 

10 

2.8 

MSA 

1.4 

1.2 

1.2 

28 

10 

2.8 

ITERATIVE  SCA 

WITHOUT  SUBARRAYS 

2.9 

1.5 

1.9 

30 

11 

2.7 

WITH  SUBARRAYS 

2. 8-3. 2 

1.5-2 .2 

1.7-1. 4 

30 

11 

2.7 

NONITERATIVE  SCA 

WITHOUT  SUBARRAYS 

52-250 

6-23 

9-11 

79-277 

15-32 

7-9 

WITH  SUBARRAYS 

11-43 

1.9-2. 5 

6-17 

38-70 

11-12 

3.5 

TABLE  5.1  RELATIVE  CPU  TIME.  DSA  IS  FASTEST  FOR  BOTH  1-  AND  2-D  IMAGES.  NON¬ 
ITERATIVE  SCA  IS  SLOWEST.  ITERATIVE  SCA  WITH  OR  WITHOUT  SUBARRAYS 
IS  HIGHLY  EFFICIENT. 

a.  Limerick  Nuclear  Power  Plant,  83-m  array. 

b.  Boeing  727,  ISAR  Image.  Target  is  in  the  clear. 

c.  Ratio  of  a  to  b. 
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6 .  SUMMARY 

It  is  shown  theoretically  and  experimentally  corroborated 
that  highly  distorted  phased  arrays  can  be  self-calibrated  under  a 
very  wide  range  of  conditions.  The  dominant  scatterer  algorithm 
(DSA) ,  which  is  the  basic  procedure  in  the  first  class,  can  be 
made  to  self-calibrate  almost  always  on  man-made  structures.  A 
modification,  the  multiple  scatterer  algorithm  (MSA),  extends  its 
usefulness  under  conditions  of  poor  target  to  clutter  ratio.  The 
scene-free  ABF  algorithms  are  superior  when  the  source  distribu¬ 
tion  is  statistically  homogeneous,  e.g.,  sea  clutter.  Although  the 
scene-free  techniques  were  designed  primarily  for  statistically 
homogeneous  backscattering,  the  experimental  evidence  indicates 
that  they  perform  exceedingly  well  under  those  conditions  in  which 
the  dominant  scatterer  class  of  algorithms  work.  Thus  the  scene- 
free  algorithms  have  been  found  to  be  more  generally  applicable 
than  the  dominant  scatterer  algorithms .  They  do,  however,  require 
greater  computational  strength.  However,  because  of  the  continuing 
improvement  over  the  years  in  speed,  capacity  and  cost  of  realtime 
signal  processing  devices,  it  is  concluded  that  the  design  of 
future  radar  systems  requiring  very  large  antennas  should  be  based 
upon  the  scene-free  algorithms. 

The  reader  is  referred  to  pages  103,  158,  182  and  183  of  the 
Appendix  for  more  detailed  summary  observations. 
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